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Abstract. We give a new heuristic for all of the main terms in the integral moments of 
various families of primitive L-functions. The results agree with previous conjectures for 
the leading order terms. Our conjectures also have an almost identical form to exact ex- 
pressions for the corresponding moments of the characteristic polynomials of either unitary, 
orthogonal, or symplectic matrices, where the moments are defined by the appropriate group 
averages. This lends support to the idea that arithmetical L-functions have a spectral inter- 
pretation, and that their value distributions can be modeled using Random Matrix Theory. 
Numerical examples show good agreement with our conjectures. 
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1. Introduction and Statement of Results 

Random Matrix Theory (RMT) has recently become a fundamental tool for understand- 
ing L-functions. Montgomery |Monj showed that the two-point correlations between the 
non-trivial zeros of the Riemann ^-function, on the scale of the mean zero spacing, are sim- 
ilar to the corresponding correlations between the eigenvalues of random unitary matrices 
in the limit of large matrix size Mehj and conjectured that these correlations are, in fact, 



identical to each other. There is extensive numerical evidence jOdlj in support of this conjec- 
ture. Rudnick and Sarnak |RSj extended Montgomery's analysis to all n-point correlations, 
and to the zeros of other principal L-functions. Katz and Sarnak |KSaj introduced the idea 
of studying zero distributions within families of L-functions (see also |OS[ IRub] ) and have 
conjectured that these coincide with the eigenvalue distributions of the classical compact 
groups. In this context symmetries of an L-function family determine the associated classi- 
cal group. We shall here be concerned with the distribution of values taken by L-functions, 
either individually (i.e. along the appropriate critical line), or with respect to averages over 
families. Specifically, we shall calculate the integral moments of these distributions. 

Keating and Snaith [KSl j suggested that the value distribution of the Riemann (- function 
(or any other principal L-function) on its critical line is related to that of the characteristic 
polynomials of random unitary matrices. This led them to a general conjecture for the 
leading-order asymptotics of the moments of this distribution in the limit of large averaging 
range. Their conjecture agrees with a result of Hardy and Littlewood |HLj for the second 
moment and a result of Ingham jl] for the fourth moment (see, for example [Tj). It also agrees 
with conjectures, based on number-theoretical calculations, of Conrey and Ghosh |CG2j and 
Conrey and Gonek |GGoj for the sixth and eighth moments. General conjectures for the 
leading-order asymptotics of the moments of L-functions within families, based on random- 
matrix calculations for the characteristic polynomials of matrices from the orthogonal and 
unitary-symplectic groups, were developed by Conrey and Farmer |GFj and Keating and 
Snaith |KS2j . These are also in agreement with what is known, and with previous conjectures. 

Our purpose here is, for the integral moments of a family of primitive L-functions, to go 
beyond the leading order asymptotics previously investigated: we give conjectures for the 
full main terms. We propose a refined definition of "conductor" of an L-function, which 
to leading order is the (logarithm of) the "usual" conductor. We find that often, but not 
always, the mean values can be expressed as polynomials in the conductor. Importantly, our 
conjectures show a striking formal similarity with analogous expressions for the character- 
istic polynomials of random matrices. This provides a strong measure of the depth of the 
connection between L-functions and RMT. We also perform numerical calculations which 
show very good agreement with our conjectures. Non-primitive families can also be handled 
by our methods, but we do not treat those here. 

The conjectures we develop here can also be obtained by techniques of multiple Dirichlet 
series, as described by Diaconu, Goldfeld, and Hoffstein |DGHj . In their formulation, one 
considers Dirichlet series in several complex variables. The mean values we conjecture would 
then follow from a plausible conjecture about the polar divisors of the function. An interest- 
ing feature of their approach is that for higher moments it seems to predict lower order terms 
of the form cT A with | < A < 1, while in this paper we conjecture that our main terms are 

valid with an error of size 0{T^ +£ ). The cubic moment of quadratic Dirichlet L-functions is 
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a specific case for which there is a conjectured lower order term fZ. which possibly could be 
tested numerically. 

There are many theorems dealing with moments of L-functions in particular families. 
The technique to prove these theorems usually involves invoking an approximate functional 
equation and averaging the coefficients of the L-function over the family. The averaging 
process behaves like a harmonic detection device. This harmonic detector usually presents 
itself as a formula with a relatively simple part and a somewhat more complicated part that 
is smaller in the first approximation. In the theorems in the literature it is often the case 
that the simple part of the harmonic detector is sufficiently good to determine the first or 
second moment of the family. The terms involved here are usually called the "diagonal" 
terms. But invariably the more complicated version is needed to determine the asymptotics 
of the third or fourth moments; in these situations one has gone "beyond the diagonal." In 
at least one situation (fourth moment of cusp form L-functions) it has been necessary to 
identify three separate stages of more subtle harmonic detection: the first featuring diagonal 
term contributions and the second and third featuring contributions to the main terms by 
two different types of off-diagonal terms. We believe that as one steps up the moments of a 
family then at every one or two steps a new type of off-diagonal contribution will emerge. 
The whole process is poorly understood; we only have glimpses of a mechanism but no clear 
idea of how or why it works. 

It is remarkable that all of these complicated harmonic detection devices ultimately lead 
to very simple answers, as detailed in this paper. It is also remarkable that there are only 
three or four different types of symmetries; families with the same symmetry type often 
have different harmonic detectors, with different wrinkles at each new stage of off-diagonal, 
but somehow lead to answers which are structurally the same. It would be worthwhile to 
understand how this works. 

Finally, we comment that the recipe we develop in this paper only uses the simplest 
diagonal harmonic detectors. Our formulas are expressed as combinatorial sums arising only 
from diagonal terms. We are well aware of the off-diagonal pieces, and we do not understand 
how they cancel and combine. What we do understand and what we are presenting here 
is a conjecture for the final simple answer that should emerge after all of the complicated 
cancellations between the increasingly subtle off-diagonal terms are taken into account. The 
reader needs to be aware of this to understand the goals and contents of this paper. 

The paper is organized as follows. In the remainder of this section we give a detailed 
comparison between L-functions and characteristic polynomials of unitary matrices, summa- 
rize our previous work on the leading terms in the mean values of L-functions, and describe 
the more general moments considered in this paper. This allows us to state our main results 
and conjectures, which are given in Section 11.51 We then give a detailed comparison with 
known results for the Riemann (^-function. 

In Section|21we give a detailed derivation of our conjectures in the case of moments on the 
critical line of a single L-function. We first write the conjecture in terms of a function defined 
by an infinite sum, and then write it as an Euler product and identify the leading-order poles. 
The local factors are also written in a concise form which is more suitable for computation. 
Both the L-function and random matrix calculations lead to expressions involving a sum 
over a set of partitions. These sums can be written in a concise form involving contour 
integrals, as described in Section 12.51 We also show that the original results of Keating and 
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Snaith [KSll lK^2] for the leading order term can be re-derived from the present work. In 
addition, we express the arithmetic factor in the moments of the Riemann zeta-function in 
an explicit form. 

In Section El we describe a particular notion of a family of L-functions which can be used 
to give a unified treatment of all of the mean values we have considered. These families 
are central to our method of conjecturing mean values and we give a detailed description 
of the method in Section HJ As explicit examples we give the details of the calculations for 
L-function families with Unitary, Symplectic, and Orthogonal symmetry. 

In Section|S]we give numerical approximations for the coefficients in our conjectured mean 
values. We then report on numerical calculations of representative cases of the conjectures. 
Good agreement is found. 

The calculations of the random matrix averages, which are based in part on [BF] and 
BH , are complicated but elementary. Those results have been presented in [CFKRS . In 
subsequent papers we will also present a fuller discussion of the terms which appear in our 
conjectures, give some more general conjectures, and describe the algorithms behind our 
numerical calculations. 

The authors thank P. Forrester, R. Heath-Brown, C. Hughes, N. Katz, P. Michel, and 
P. Sarnak for many helpful discussions. 

1.1. Properties of L-functions. We present the definition and key properties of L-functions 
These properties are familiar, but a summary will be useful in our discussion of mean values 
and for the comparison with the characteristic polynomials of random matrices. 

The definition of an L-function which we give below is a slight modification of what has 
come to be called the "Selberg class" ICG31 IMurj of Dirichlet series. Let s = o + it with 
a and t real. An L-function is a Dirichlet series 

oo 
n=l 

with a n <C £ n e for every e > 0, which has three additional properties. 

Analytic continuation: L(s) continues to a meromorphic function of finite order with at 
most finitely many poles, and all poles are located on the a = 1 line. 

Functional equation: There is a number e with |e| = 1, and a function Jl( s ) of the form 

w 

(1.1.2) 7L (s)=P(s)Q s l[T(w,s + fi J ), 

i=i 

where Q > 0, Wj > 0, K/i,- > 0, and P is a polynomial whose only zeros in a > are at the 
poles of L(s), such that 

(1.1.3) t L (s):= lL {s)L{s) 
is entire, and 

(1-1-4) Us)=eUl-s), 

where £l(s) = £l(s) and s denotes the complex conjugate of s . 
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The number 2 Y^j=i w j * s called the degree of the L-function, and this is conjectured to 
be an integer. It is conjectured furthermore that each Wj can be taken to equal |, so w 
equals the degree of the L-function. 

For the calculations we do in this paper, it is convenient to write the functional equation 
in asymmetric form: 

(1.1.5) L(s)=eX L (s)L(l- 8 ), 
7l(1 — s ) 

where X L (s) = — — . Also we define the "Z-function" associated to an L-function: 

(1.1.6) Z L (s):=e^XP(s)L(s), 
which satisfies the functional equation 

(1.1.7) Z L (s) = Z2(l-s). 

Note that here we define as a function of a complex variable, which is slightly different 
from the standard notation. Note also that Zl{\ + it) is real when t is real, -Xx(f) = 1' anc ^ 
\X L {\ +it)\ = 1 if t is real. 
Euler product: For o > 1 we have 

(1.1.8) L(s) = l[L p (l/p s ), 

v 

where the product is over the primes p, and 

(1-1-9) W) = E? = exptel , 

where b n <C n 6 with 9 < ^. 

Note that L(s) = 1 is the only constant L-function, the set of L-functions is closed 
under products, and if L(s) is an L-function then so is L(s + iy) for any real y. An L- 
function is called primitive if it cannot be written as a nontrivial product of L-functions, 
and it can be shown, assuming Selberg's orthonormality conjectures, that any L-function 
has a unique representation as a product of primitive L-functions. See |CG3j . It is believed 
that L-functions only arise from arithmetic objects, such as characters II )av| . automorphic 
forms |Iwlt IIw2j . and automorphic representations |BCt IBuj . Very little is known about 
L-functions beyond those cases which have been shown to be arithmetic. 

There are several interesting consequences of the above properties, and many conjectures 
which have been established in a few (or no) cases. We highlight some additional properties 
of L-functions and then discuss their random matrix analogues. 

Location of zeros: Since £l(s) is entire, L(s) must vanish at the poles of the T-functions 
in the 7^ factor. These are known as the trivial zeros of the L-function. By the functional 
equation and the Euler product, the only other possible zeros of L(s) lie in the critical strip 
< a < 1. By the argument principle, the number of nontrivial zeros with < t < T is 
asymptotically (W/n)T\ogT, where W = Yl w j- The Riemann Hypothesis for L(s) asserts 
that the nontrivial zeros of L(s) lie on the critical line o = |. The much weaker (but still 
deep) assertion that L(s) 7^ on o = 1 has been proven for arithmetic L-functions |JSj . 
which can be viewed as a generalization of the prime number theorem. 
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Average spacing of zeros: By the zero counting result described above, the average gap 
between consecutive zeros of L(s) with imaginary part around T is 7r/(WlogT). 

Zeros of derivatives: If the Riemann Hypothesis is true then all zeros of the derivative £'(s) 
lie on the critical line, while all zeros of C{s) ^ e to the right of the critical line |LM| . 

Critical values: The value L(f) * s cane d the critical value of the L-function. The significance 
of s = \ is that it is the symmetry point of the functional equation. The mean values we study 
in this paper are averages of (powers of) critical values of L-functions, where the average is 
taken over a "family" of L-functions. Examples of families and their corresponding mean 
values are given in Section IQ1 

Note. If the set {/jlj} is stable under complex conjugation and the a n are real, then e is 
commonly called the sign of the functional equation. If the sign is —1 then L(s) has an odd 
order zero at s = |; more generally, if the sign is not 1 then L(|) = 0. When L(|) vanishes, 
it is common to use the term 'critical value' for the first nonzero derivative L^(|), but in 
this paper we use 'critical value' to mean 'value at the critical point.' 

Log conductor: We measure the "size" of an L-function by its log conductor, defined as 
c(L) = cond(L) = \X' L (^)\. The conductor of an L-function has a conventional meaning 
in many contexts, and the log conductor is a simple function of the (logarithm of the) 
usual conductor. Other authors use similar names, such as "analytic conductor" , for similar 
quantities. By the argument principle, the density of zeros near the critical point is 2irc(L)~ 1 . 

Approximate functional equation: A standard tool for studying analytic properties of L- 
functions is an approximate functional equation for L(s), which expresses the L-function as a 
sum of two Dirichlet series involving the Dirichlet coefficients of L multiplied by a smoothing 
function. See for example |1K| 5.3]. For the purposes of the heuristics we develop, we use a 
sharp cutoff and don't concern ourselves with the remainder, 

(1.1.10) L(s) = — + eX L (s) — - — h remainder. 

Here the product xy depends on parameters in the functional equation. The name comes 
from the fact that the right side looks like L(s) if x is large, and like eXl(s)L(1 — s) if x is 
small, which suggests the asymmetric form of the functional equation. 

The approximate functional equation is the starting point of our approach to conjecturing 
the moments of L-functions. This is described in Sections 12. II and 14. II 

1.2. Properties of characteristic polynomials. With the exception of the Euler product, 
all of the properties of L-functions have a natural analogue in the characteristic polynomials 
of unitary matrices. We note each property in turn. 

Let 

N 

(1.2.1) A(s) = A A (s) = det(J - A*s) = JJ (l - se~ Wn ) 

n=l 

denote the characteristic polynomial of an N x N matrix A. Throughout the paper we 
assume that A is unitary (i.e. A* A = I where A* is the Hermitian conjugate of A), so the 
eigenvalues of A lie on the unit circle and can be denoted by e tdn . 
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Note: in our previous paper CFKRS we used a different definition of the characteristic 
polynomial. 

We can express A(s) in expanded form: 

N 

(1.2.2) A(s) = J> n s", 

n=0 

which corresponds to the Dirichlet series representation for L-functions. 
Analytic continuation: Since A(s) is a polynomial, it is an entire function. 
Functional equation: Since A is unitary, we have 

(1.2.3) A A (s) = {-l) N detA* s N det(J - As' 1 ), 
and so, writing 

(1.2.4) detA = e itt> 
(where unitarity implies that <fi £ R), we have 

A A (s) = (-l) N detA*s^A A ^ s ) 

(1.2.5) = (-l) 7V e-VA A (i). 

This plays the same role for A(s) as the functional equation for L-functions: it represents a 
symmetry with respect to the unit circle (s = re ia — > j = ^e~ ia ). 

Also let 

(1.2.6) Z A {s) = ((-l) N e-^y l > s- N / 2 A A (s), 

in direct analogy to (jl.l.6j) . the sign of the functional equation e being identified with 
(—l) N e~ 1 ^ = (— 1)^ det A*. The functional equation becomes 

(1.2.7) Z A (s)=Z A (\). 

Note that this implies Z is real on the unit circle, and in analogy to the Xl factor from an 
L-function, the factor s~ N ^ 2 equals 1 at the critical point s — 1, and has absolute value 1 on 
the unit circle. 

Location of zeros: Since A is unitary, its eigenvalues all have modulus 1, so the zeros of A(s) 
lie on the unit circle (i.e. the Riemann Hypothesis is true). The unit circle is the "critical 
line" for A(s). 

Average spacing of zeros: Since the N x N matrix A has eigenvalues on the unit circle, 
the average spacing between zeros of A A (s) is 2n/N. 

Zeros of derivatives: Since the zeros of A(s) lie on the unit circle, the zeros of the derivative 
A' lie inside the unit circle. This follows from the general fact that the zeros of the derivative 
of a polynomial lie in the convex hull of the zeros of the polynomial. 

Critical values: The critical point for A(s) is the symmetry point of the functional equation 
s — 1 — e*'°, and A(l) is the critical value. 

Conductor: In analogy with the case of L-functions, we define the conductor of A to be 
the (absolute value of the) derivative of the factor in the asymmetric form of the functional 
equation, evaluated at the critical point s = 1. That is, the conductor of A is N. Also in 
analogy to the case of L-functions, the density of zeros on the unit circle is 2tt/N . 
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When modeling a family of L-functions, we choose N so that L and A have the same 
conductor. Equivalently, L and A have the same density of zeros near the critical point. 

Approximate functional equation: Substituting the polynomial (jl.2.2|) into the functional 
equation (jl.2.5|) . we have 

N N 

(1.2.8) ^a nS ™ = (-l^e-^a^- 

n=0 n=0 

and so 

(1.2.9) a n = (-l) N e'^aW^. 
Hence, when A" is odd, we have 



N-l JV-1 



;i.2.10) A(s) = a ™ sm + (-l)^e-^ J2 



a n S , 

m=0 n=0 

which corresponds to the approximate functional equation for L-functions. When A" is even, 
there is an additional term: on_s~ . 

2 

Although we use the approximate functional equation in our calculations for L-functions, 
in our previous paper [CFKRSJ we use other methods for the characteristic polynomials. In 
principle, it would be possible to use the approximate functional equation and compute 
averages of products of the coefficients a n . Such a calculation would, presumably, mirror 
that for the L-functions. This would appear to be more cumbersome than the approach 
taken in |CFKRSj . but might merit further investigation. 

The above discussion applies to any unitary matrix. We also consider matrices which, in 
addition to being unitary, are either symplectic or orthogonal. We use these three ensembles 
of matrices to model families of L-functions. While the notion of "family of L-functions" 
has not yet been made precise, we give several natural examples in the next section. 

Associated to each family is a "symmetry type" which identifies the matrix ensemble 
that will be used to model the family. This correspondence is most easily seen in terms of 
the sign of the functional equation, which is analogous to the determinant of the matrix. If A 
is unitary symplectic, then det A = 1 (i.e. = 0), and if A is orthogonal, then det A = ±1. 
Correspondingly, the functional equations for L-functions with unitary symmetry involve a 
(generally complex) phase factor, whereas for L-functions with symplectic symmetry this 
phase factor is unity, and in the case of orthogonal symmetry it is either +1 or —1. 

While the sign of the functional equation can sometimes suggest the symmetry type of 
the family, in general it requires a calculation to determine the symmetry type. One possible 
calculation is to determine the moments of the family near the critical point, as described 
in this paper. Comparison with the corresponding random matrix average can then be used 
to determine the symmetry type. Another possibility is to determine the density of the 
low-lying zeros of the family. 

1.3. Example families and moments of L-functions. We now give examples of families 
of primitive L-functions and describe the associated mean values. The families we consider 
here are of a special form, which is described in Sectional In preparation for the comparison 
with random matrices in the next section, we will classify the example families according to 
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their symmetry type: Unitary, Orthogonal, and Symplectic. For the Orthogonal symmetry 
type we recognize three cases: SO, 0~, and O, corresponding respectively to Orthogonal 
families in which the functional equation has e = 1, or e = —1, or e = ±1 equally often. 
Note that each family is a partially ordered set, and the order is determined by a quantity 
called the "conductor" of the L-function. The mean values given below are conjectural for 
all but a few small values of k. For a general discussion of these mean values and some more 
examples, see |CFj . 

Unitary examples: 

1) {L(s + iy) | y > 0}, ordered by y, where L(s) is any primitive L-function. These are the 
only known continuous families of L-functions (Sarnak's rigidity conjecture). 

2) {L(s,x) | Q a positive integer, x a primitive character modg} ordered by q. 
An example conjectured mean value for integer k is: 

(1.3.1) f |C(| + it)\ 2k dt = TV k {\ogT) + 0(T^ +£ ), 

J o 

for some polynomial Vk of degree k 2 with leading coefficient gk(ik/k 2 \, where 



Tim + k) 

a k = I I I 1 - t I } A J P 



n(i-o e 



m=0 



m\T{k) 



-m 



(fc-i) 2 ^ fk - 1 



3=0 



J 



and 

(1.3.3) 9k = k 2 \\l 



k-l 



(The placement of fc 2 ! is to ensure that is an integer jCFj .) The above conjecture has 
been proven for k = 1, 2. See [HD HI El IKS IH-B11 O lMot] . When k — 2 our conjectured 
error term of 0(T2 +e ) has only been obtained in the case of a smooth weight function jlvlj . 

Symplectic examples: 

3) {L(s,Xd) | d a fundamental discriminant, Xd(n) = (-)} ordered by \d\. 

4) {L(s, sym 2 f) \ f G S , fc(ro(l)), k a positive even integer}, ordered by k. 
An example conjectured mean value is: 

(1-3.4) L (i ^) fc = ^ D Qk(logD) + 0(D^ £ ), 

|d|<D 

where ^ s over fundamental discriminants, Xd(n) = (-) is the Kronecker symbol, and 
the sum is over all real, primitive Dirichlet characters of conductor up to I*. Here Qk is a 
polynomial of degree kik + l)/2, with leading coefficient gkdk/(k(k + l)/2)!, where 

(1-3.5) a ^ll 1+ i I 2 + ; 
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and 

k .. 



1.3.6) ^ = (fc(fc + l)/2)!jl 



The main term of this conjecture has been proven for k = 1, 2, 3, and the case of k = 4 is 
almost within reach of current methods. See |GV| []J Eoj . 

Orthogonal examples: 

5) {L(s, f)\fe Sk(T (N)), N fixed, k a positive even integer}, ordered by k. 

6) {L(s, /) | / G Sk(T (N)), k fixed, N a positive integer}, ordered by N. 
An example conjectured mean value is: 

(1-3.7) J2 L f {\j) k = q -n k {\o gq ) + 0{q^% 

where #2(9) is the collection of Hecke newforms of weight 2 and squarefree level g. Here TZk 
is a polynomial of degree k(k — l)/2, with leading coefficient gy.a^/{k{k — l)/2)!, where 



1.3.8) a* = n(l"j) 2 ^/ 

_ L V 



and 



gi0 g— «6 



fc-1 



1.3.9) ff* = 2*- 1 (*(A: - 1)/2)I JJ 



11 W 

The main term of this conjecture has been proven for k = 1, 2, 3, 4, in the case that q is 
prime. See jDj IDFH IKMVj . Also see Ivic |Iv2j for the analogous mean values for Maass 
forms. 

The above examples are merely meant to give a flavor of the types of families which are 
of current interest. 

The above cases, and their random matrix analogues, have been extensively discussed 
from the perspective of the leading terms in the asymptotic expansions. See [OF] IKS1[ IKS2] . 
In the present paper we extend that work to include all of the terms in the above mean 
values (ie., all coefficients in the conjectured polynomials), which we recover from a more 
general mean value involving a product of L-functions whose arguments are free parameters. 
In the next two sections we describe these more general mean values, discuss their random 
matrix analogues, and then state our results and conjectures. 



1.4. Shifted moments. A key point in this paper is that the structure of mean values of 
L-functions is more clearly revealed if one considers the average of a product of L-functions, 
where each L-function is evaluated at a location slightly shifted from the critical point. The 
example mean values given in the previous section can be obtained by allowing the shifts to 
tend to zero. 
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Let a = (ai, . . . ,«2fc), where throughout the paper we assume |3tay| < |, and suppose 
that g(t) is a suitable weight function. The mean values we consider are 

/oo 
Z L {\ + a 1 + it)--- Z L (\ + a 2k + it)g(t) dt, 
-oo 

and, with a = (ax, . . . , «&), 

(1.4.2) S h (T, a,g) = J2 z d\ + «i) • • • Z L (\ + a fe )<?(c(L)). 

In the first case it is assumed that L(s) is a primitive L-function, and in the second T is 
family of primitive L-functions partially ordered by log conductor c(L). 

We refer to g as a "suitable" weight function, but we leave that term undefined. An ex- 
ample of a suitable weight function is g(x) = f(x/T), where / is real, nonnegative, bounded, 
and integrable on the positive real line. 

The random matrix analogs of the above expressions are 



(1.4.3) Jk(U(N), a) = / Z A (e-^) ■ ■ ■ Z A ( e - a ^) dA, 

Ju(N) 

where a = (ax, • • • , «2fc) and the average is over Haar measure on U(N), and 

(1.4.4) J k (G(N),a)= [ Z A (e-°*)---Z A (e- a *)dA, 

Jg(n) 

where G(N) is USp(2N), 0~(2N), or SO(2N) and a = (a Xj . . .,a k ). 0~(2N) is defined as 
the collection of orthogonal 2N x 2N matrices with determinant negative one. Haar measure 
on USp(2N) and 0(2N) determines the weighting for the averages. 

In the next section we compare our conjectures for the L-function mean values with 
exact formulae for the random matrix averages. 

1.5. Main results and example conjectures. We state our main results and conjectures 
here. We give example conjectures for the full main term in shifted mean values of number 
theoretic interest; these examples illustrate our methods and cover the three symmetry types 
of families of L-functions. We also give a corresponding theorem about the random matrix 
analogue of these mean values for each of the three compact matrix ensembles. 

We present our results in pairs: a conjecture for an L-function mean value, followed by 
a theorem, quoted from |CFKRSj . for the corresponding average of characteristic polyno- 



mials. For each pair the parts of each formula match according to the scheme described in 
Section rOl In particular, the scaling of the large parameter is determined by equating log 
conductors. In the random matrix formula the integrand contains a term (1 — gi^™-^)- 1 ; 
which has a simple pole at zi = ±z m . In the L-function formula this corresponds to the 
term containing all of the arithmetic information, which is of the form £(1 + Zi =F z j) times 
an Euler product, and so also has a simple pole at Z; = iZj. 

The formulae are written in terms of contour integrals and involve the Vandermonde: 
(1.5.1) A(zx,...,z m )= H 

l<i<j<m 

We also set e(z) = e 2mz . 
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Conjecture 1.5.1. Suppose g(t) is a suitable weight function. Then 

/oo /»oo 
\C(i + it)\ u g{t)dt= / P k (log^)(l + 0(r^))g(t)dt, 
-OO J —OO 

where is the polynomial of degree k 2 given by the 2k- fold residue 
(1.5.3) 



fc! 2 (27Ti) 2fc / / 2fe - _ i ..._ 2fe , 

where one integrates over small circles about Z{ = 0, with 

k k 

(1.5.4) G(zi, z 2k ) = A k (zi, . . . , z 2k ) Yi II C(l + - zjfc+i), 

i=l 3=1 

and Afc is t/ie Euler product 

<-) *w - rmn - ?*M jf n (> - $Y - 

More generally, 

/oo 
P fe (log£,a)(l + 0(H +£ ))a(t) 
-oo 

w/iere 
(1.5.7) 

p / \ (~l) fc 1 j" / G(z 1} z 2fc ) A 2 (zi, . . . , z 2k ) £ £ » Zj _ Zfc+ . , , 
P k (x, a) = fc|2 fc y> • • • j — e2 ^ » <+> d Zl ... dz 2k , 

nn&-«o 

3=1 i=l 

TOt/j t/ie pai/i of integration being small circles surrounding the poles aj. 

A general version of the above conjecture is given in Conjecture 12.5.41 
Theorem 1.5.2. In the notation of Section \l-4\ we have 

k-1 / .j k \ 

(1.5.8) J k (U(N), 0) = J] \^-^ J[{N + i + j)j . 

More generally, with 

k k 

G(z 1 ,...,z 2k ) = lll[(l-e- z > +z ^r 1 , 

i=l 3=1 

we have 

i=l j = l 
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Comments on the formulae: 

(1) Let on — > in the second part of the Conjecture to obtain the first part of the 
Conjecture. 

(2) The structures of Jk(U (N), a) and Pk(x, a) are identical in that the functions G(zi, . . . . 
have simple poles at Zi = Zk+j. 

(3) The local factors of Ak(a) are polynomials in p~ x and p~ ai , i = 1, . . . , k, as seen 
from Theorem 12.6.21 Since A^a) comes from a symmetric expression, it is also a 
polynomial in i — 1, . . . ,k. This is discussed in Section l2~o1 Note also that 

in (jl.3.2j) equals Ak(0, ■ ■ ■ , 0), as shown in Section l2~71 

(4) That Pk{x) is actually a polynomial of degree k 2 can be seen by considering the order 
of the pole at z 3 ■ = 0. We wish to extract from the numerator of the integrand, the 
coefficient of \\ z i k ~ ] \ a polynomial of degree 2k(2k — 1). The Vandermonde determi- 
nant squared is a homogeneous polynomial of degree 2k(2k — 1). However, the poles 
coming from the £(1 + Zi — z^+j) cancel k 2 of the Vandermonde factors. This requires 
us, in computing the residue, to take, in the Taylor expansion of exp(| Y^i z j ~ z k+j), 
terms up to degree k 2 . 

(5) The fact that Pfc(log^r) is a polynomial in log j- of degree k 2 corresponds nicely 
to the formula for Jk(U(N),0), which is a polynomial of degree k 2 in N. Equating 
the density of the Riemann zeros at height t with the density of the random matrix 
eigenvalues suggests the familiar equivalence N = log^r. In this paper we view this 
as equating conductors. 

(6) The leading term of Pk{x) coincides with the leading term conjectured by Keating 
and Snaith (see Section l2~7j) . The full polynomial Pk{x) agrees, when k — 1 and 
k — 2, with known theorems (see Sections II .61 and II. 7|) . 

(7) We can recover the polynomial Vk in (jl.3.1j) from P^ by taking g(t) = X[o,T]{t) in the 
conjecture. 

(8) The multiple integrals in the Theorem and the Conjecture can be written as combi- 
natorial sums. See Section El where a detailed derivation of our conjecture is given. 

(9) Our conjecture concerning the order of the error term is based on our numerical 
calculations (see Section |3j) and examination of examples in the literature. 

Conjecture 1.5.3. Suppose g(u) is a suitable weight function with support in either (0, oo) 
or (—oo,0) ; and let Xd(s) = \d\ z~ s X(s, a) where a = if d > and a = 1 if d < 0, and 



'1.5.10) X(s, a) = n s - l 2 \ 2 J- 



That is, Xd( s ) is the factor in the functional equation L(s, Xd) — £dXd{s)L(l — s, Xd)- Sum- 
ming over fundamental discriminants d we have 



(1.5.11) Yl L (iXd) k g(\d\) = Y. Qk(log\d\)(l + 0(\d\-^))g(\d\) 
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where Q k is the polynomial of degree k(k + l)/2 given by the k-fold residue 
(1.5.12) 

^.^g^g i;, j o {zi , ) A W ,...,^ ef£ , i , fei dzt 

II 

i=i 

u>/iere 



(1.5.13) G(z ll ...,2 fe )=4(2 1 ,...,^)J]l(i + z i ,a)-5 J] C(l + ^ + 

and Afc zs i/ie Euler product, absolutely convergent for < |, defined by 



z j), 



A k ( Z1 ,..., Zk ) = n n (i-^-) 



P l<i<jr'<fc 

More generally, if ' T is the family of real primitive Dirichlet L-functions then 



(1.5.15) S k (F,a,g) = J2 Q k(\og \d\, a) (I + 0(\d\-^))g(\d\), 

d 

in which 

(l\k(k-l)/2 2 k i 

Q k (x,a) = — ; 



k\ (2m) k 

r i G(z u z k )A(zl zp 2 n k 3=1 z A | ^fc 

(1.5.16) x &>■■■&> — e2^ =1 J dzi . . . cfefc 

i=i i=i 

where the path of integration encloses the ±a 7 s. 
Theorem 1.5.4. In the notation of Section \l.J\ we have 

(1.5.17) J k (USp(2N),0)= (^ k(k+1)/2 Uj^) II + 

More generally, with 

G(zi, ...,z k ) = \ \ (1 - e" 

i<i<ji<fc 

we nave 

J i(f / S p(27V),a) = L_i_ _ 

/ 1K1Q x / / G ( z l> ■ ■ ■ , z k)&(z 2 , • • • , z k) 2 Il k j=i z j N yk , 2 . , , 

(1.5.18) x (p ■ ■ ■ (p — - e ^^i^dzx ■ ■ ■ dz, 



l<i<j<k 



YlY[( z j - ®i)( z j + ®i) 
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where the contours of integration enclose the ±a 's. 
Comments: 

(1) When comparing Theorem 1 1 . 5 . 41 wit h Conjecture II. 5. 3[ equating log conductors (i.e., 
the density of zeros) gives the equivalence 

(1.5.19) 2N = cond(d) := log(|d|/vr) + (r'/r)(± + a). 

The conductor we use here should be contrasted with the "usual" conductor asso- 
ciated with Dirichlet L-functions: log(|d|/7r) — log(2). We believe this difference is 
significant, so we discuss it briefly. 

The following manipulations show that our conductor arises naturaly. In the 
derivation of the conjecture, one encounters the function Xd(~ + z)~^, which can 
we rewritten in several ways: 

X d {\ + z)-* = e^ d - z X(\ + z,a)-% 

(1.5.20) = e^ cond{dyz g{z), 

where Q{z) = 1 + 0(z 3 ). In the statement of the conjecture we used the first line 
of (jl.5.2()|) . incorporating the product over X(^ + z, a)~5 into the factor G(zi, . . . , Zk). 
If we chose instead to use the second line of fll.5.2()j) . then the conjecture would be 
written as a sum over Qk{cond(d)). One would still have that Qk is a polynomial of 
degree fc(fc+l)/2. Since Q(z) = l + 0(z 3 ), the first 3 leading terms in that polynomial 
would not explicitly depend on the factor X^ from the functional equation, although 
the lower degree terms would. This phenomenon does not occur for moments of 
L-functions in t-aspect. 

(2) Qk in (jl.3.4|) can be recovered from Qk above by taking g(\d\) = X[o,D](\d\), and using 
the estimate ^ -D<d<o ^ = 3D/tt 2 + 0(D^ +6 ); the same estimate holds for positive d. 

(3) A heuristic derivation of Conjecture 11.5.31 is given in Section l4~4l 

(4) The leading term of Qk coincides with the leading term conjectured by Keating and 
Snaith |KS2j . The calculation is analogous to the one given in Section 12.71 

Conjecture 1.5.5. Suppose q is squarefree, let H n (q) be the set of new forms in S n {T (qj), 
and let 



1.5.21) X n Js)={^) T 3 f + f ) 



Air 2 J r(s-i + § 

be the factor in the functional equation Lf(s) = e n ^X n ^ q (s)Lf(l — s) for the L-functions 
associated to f e H n (q). Then 

(1.5.22) Yl L f(l) h (/J)" 1 = E Rk(n, q )(fj)- 1 (l + 0(nq)-^) 

feH n ( q ) feH n (q) 

as nq — > oo, where Rk(n,q) is given by the k-fold residue 
(1.5.23) 

Rk ~ ifej (27ri)* J "J H k 1 z 2k - 1 ll A ^lf+^J *dz x ...dz k 
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where 



Zi + Zj) 



(1.5.24) G(z 1 ,...,z k ) = A k (z 1 ,...,z k ) J] C(l + 

l<i<j<k 

and A k is the Euler product which is absolutely convergent for \$lzj\ < \, with j = 1, . . . , k, 
defined by 



. . . , z k ) = n n (i-^Tij) 



p\q l<i<j<k 

(^.25) x 2 -j o i.'»n v - w,/ p 7 

To state the more general version of Conjecture 11.5.51 involving a sum of products of 
Lf(h + Uj), it is natural also to consider the sums over even / and odd / separately. See 
Conjectures 14.5.11 and 14.5.21 

Theorem 1.5.6. In the notation of Section \l.J\ we have 

(1.5.26) J k (SO(2N), 0) = ( 2^ +1 )/ 2 J] J J] (iV + 
More generally, with 

(1.5.27) G(z 1 ,...,« fc )= II (l-e-") _1 

l<^<m<fc 

we have 

(_l)k(k-l)/2 2 k 1 

(1.5.28) x / • • • / gfa,- t -,^,-,^n^ V ^ « idzi . . . ^ 

JJ -Oi)(«j + ai) 

i=i j=i 

and 

f_ll*(*-l)/20* 1 

J k (0-(2N),a) = i~ k 



k\ (2m] 



k 



;i. 5 .29) x # ■ . . # Cte---^,-,4) 8 I^,° v ^,^ t . . . dZk 

Y[ -OLi)(Zj +Cii) 

i=l 3=1 

Comments: 

(1) R k (n,q) does not actually depend on / e H n (q). We write (jl.5.22|) in this manner 
to stress that R k (n,q) is the expected value of Lj(|) fc . 
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(2) To compare Theorem 1 1 . 5 . 61 and Conjecture [ESS equating conductors gives the equiv- 
alence 

2N = cond(n, q) := log(g/4vr 2 ) + (T'/T)(n/2) 

(1.5.30) = \og(qn/8n 2 ) + 0(n- 1 ). 

One can express the conjectured mean value in terms of the conductor in the following 
way. In (jl.5.231) we can write 

(1.5.31) X n , q {\ + Zj )^ = e^^Gizj), 

where Q(zj) = 1 + 0(z 3 ). As in Conjecture ll.5.3| we can express Rk(n,q) as a 
polynomial in the conductor, the first 3 terms of which do not depend on the X n ^ q 
factor in the functional equation. 

(3) All of our conjectures naturally contain a factor of the form n^(| ^ z j) ~ 5 > ft J us t 
happens that in some cases A(| ± Zj) can be closely approximated by a simple 
function of the conductor. It is interesting that this same factor occurs in all of the 
random matrix moments. In that case X(s) = s~ M , where M = N or 2JV, so in the 



formula for the moments there occurs n^( e±Zj ) 2 = e 2 ^ ± 



■j 



1.6. 2nd Moment of the Riemann zeta-function. Now we consider the 2nd moment of 
the Riemann zeta-function in detail, putting our results in the context of the literature. 

Ingham's result jl] on the 2nd moment can be stated as 
(1.6.1) / C(s + a)£(l-s-0)dt= [ (C(l + a-P)+r (S - a C(l + p-a))(l + 0(r^))dt 



where s = | + it and r = r(t) = this is valid for \a\, \f3\ < |. If we let a and (3 approach 
here we obtain Ingham's theorem 

(1.6.2) / T |C(|+*t)| 2 dt= f T {\og± + 2 1 )dt + 0{T^). 
Jo Jo 

Our conjecture is compatible with these results, because, when k — 1, the function G(ai, 02) 
that appears in Conjecture 11.5. II equals £(1 + «i — a 2 ). Computing the residue, we find 

P x {x) = ± j I C(1 + Zl ~*f Z2 ~ Zl)2 e^-^d Zl dz 2 

(1.6.3) = x + 2 7 . 

The second moment with a different weighting is now given; this theorem is a slight variation 
of the theorem of Kober presented in Titchmarsh [Tj and was inspired by the numerical 
calculations described in Section I5~T1 

Theorem 1.6.1. Let 

POO 

(1.6.4) I(a,P,6)= / t(\ + it + a)({\-it-f3)e- 6t dt. 

Jo 

Then, for any r\ > 0, \a\, \/3\ < | — r] and \ arg5| < | — r] , we have 
(1.6.5) 

I(a,/3,S)= K(l + a-{3)+ (J^j C(l-a + f3)je- 5t dt + C s (a,(3) + 0(6\ogl/6) 
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uniformly in a, [3, and 5 where Cg(a,j3) <C log 1/5 uniformly in a and [3 and where 
C s (a,-a) = -2-k((2o). 

We restate the case f3 = —a as 
Corollary 1.6.2. For any fixed a with \a\ < | ; we have 



lim^ \((l + a + tt)\ 2 e- st dt - ^(l + 2ar)+Q^ ((1 - 2a) j e~ St dt 



= -n((2a) - (27r) 2a C(l - 2a)T(l - 2a) sinvra 

(1.6.6) = -27rC(2a). 

Note that 

(1.6.7) lim(C(l + x) + w- x T(l - x)C(l - x)) = (2 7 + logw). 
Thus, letting a — > 0, gives 

Corollary 1.6.3. 



1.6.8) lim QT |C(| + it)\ 2 e-» dt - J™ ^2 7 + log (J-Jj e~ St dt"} 



71". 



Remark: We discovered this Corollary after seeing the numerical results of Section 15.11 
This result also follows from a result of Hafner and I vie IHII . 



1.7. 4th Moment of the Riemann zeta- function. Now we consider the 4th moment of 
the Riemann zeta-function in detail. Our discussion here builds upon work of Atkinson |A"|. 
Heath-Brown jH-Blj . ConreyjUj, and Motohashi pot] . 

Examining Motohashi's results in detail, consider 

/oo 
C(s + Ul ) . . . C(s + u k )((l - 8 + v 1 )...£(l-s + v k )g{i) dt 
-oo 

for a function g(t) which is analytic in a horizontal strip \^s{t)\ < c and decays sufficiently 
rapidly. Motohashi obtains an exact formula for these moments for k = 1 and k = 2. We 
reformulate Motohashi's theorem (k = 2) in our context. Let 

(1.7.2) C(v) = J^X- 
v ' K ' 2 cos ft; 

and let 

(1.7.3) G s (^) = ^|. 

T(s — v) 

Then, in a notation analogous to Motohashi's, the k = 2 case of (|1.7.1|) equals 

(1.7.4) L r + L d + L c + L h , 

where L r is the (residual) main term which we are interested in here: 

POO 

(1.7.5) L r (u,v)= / W{t,u,v)g(t)dt, 



20 J. B. CONREY, D. W. FARMER, J. P. KEATING, M. O. RUBINSTEIN, N. C. SNAITH 

with 

W{t,u,v) = C{0){G s {0,0) + G 1 _ s {0,0))Z{u 1 ,u 2 ,v 1 ,v 2 ) 

+ C(lll + V 1 )(G s (u 1 ,V 1 ) + Gi- s (Ul,Vi))Z(-V 1 ,U2, -ui,v 2 ) 

+ C(ui + v 2 )(G s (u 1 ,v 2 ) + G 1 - s (u 1 ,v 2 ))Z(-v 2 ,u 2 ,v 1 , -Mi) 
+ C{u 2 + v 1 )(G s {u 2 , vi) + Gi_ s (m 2 , v i))Z(u!, -vi, -u 2 , v 2 ) 
+ C(u 2 + v 2 )(G s (u 2 ,v 2 ) + Gi- 8 (u 2 ,V2))Z(u lj -v 2 ,vi, -u 2 ) 
+ C(«! + u 2 + v 1 + v 3 )(G a (u 1 , vx)G s {u 2 , v 2 ) + G^ s {u u v 1 )G l _ s {u 2 , v 2 )) 

(1.7.6) xZ(-vi, -v 2 , -ui, -u 2 ), 

where s = \ + it, and 

1.7.7 Z{U 1 ,U 2 ,V 1 ,V 2 ) = T77— ; ; ; r • 

C(2 + Mi + u 2 + v 2 + v 2 ) 

This formula may be obtained from Motohashi's work ( |Mot] pp. 174 - 178) by a careful 
analysis of his terms together with appropriate use of the functional equation in the form 

and some trigonometric identities. 
If we use the approximation 

(i-7.9) ^^ = mr-"(i+o(i/\ s \), 

we have, using r = |t|/(27r), 

(1.7.10) C(u){G s {u) + = r- u (l + 0(l/r)) 
and 

(1.7.11) C{u + v)(G s (u)G s (v) + Gi-s(u)Gi- 8 (v)) = t' u - v (1 + 0(l/r)). 
We then have 

W{t, u, v) = (Z{u x , u 2 , vi, v 2 ) + t~ Ui - Vi Z(-Vi, u 2 , -mi, v 2 ) + r- Ul - V2 Z(-v 2 , U 2 , Vi, -Mi) 

+ T- U2 -^Z(U U -V U -U 2 , V 2 ) + T- U ^Z(U U -V 2 , Vl , ~U 2 ) 

(1.7.12) + r-^-^Zi-Vt, -v 2 , -mi, -u 2 j) (1 + 0(l/r)). 



By the formulae in Sections 12.21 and 12.51 the above agrees with the k = 2 case of Conjec- 
ture nxu 

A residue computation shows that our conjecture can be restated as 
"T pT 1 1 r r / + \ x +v 



where we integrate around small circles centered on the origin. This is in contrast to Con- 
jecture 11.5.11 which when k = 2 expresses the formula in terms of four contour integrals. It 
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may be that our formulae can be similarly simplified for all k, but we have not succeeded in 
doing so. 

2. Moments in £-aspect 

The principle behind our method of conjecturing mean values is that the Dirichlet series 
coefficients of L-functions have an approximate orthogonality relation when averaged over a 
family. These orthogonality relations are used to identify the main terms in the mean values. 

In this section we give a detailed account of the case of moments of a single primitive 
L-function. We describe the recipe for conjecturing the mean values, applying it first to 
the case of the Riemann ^-function, and then to a general primitive L-function. In the 
remainder of this section we manipulate the formulas into a more usable form, and also 
obtain a generalization of Conjecture 11.5.11 Later in Section |3] we recast our principles in a 
more general setting and consider the averages of various families of L-functions. 

2.1. The recipe. The following is our recipe for conjecturing the 2fcth moment of an L- 
function: 

(1) Start with a product of 2k shifted L-functions: 

(2.1.1) Z(s, a 2k ) = Z{\ + «!)••■ Z{\ + a 2k ) 

(here we have written the Z-function, but the examples below will show that the 
method applies to either the L- or the Z-function). 

(2) Replace each L-function with the two terms from its approximate functional equation, 
ignoring the remainder term. Multiply out the resulting expression to obtain 2 2k 
terms. 

(3) Keep the ) terms for which the product of x-factors from the functional equation 
is not rapidly oscillating. Use (J2.2.9|) to simplify the nonoscillating ^-factors. 

(4) In each of those ( k ) terms, keep only the diagonal from the sum. 

(5) Complete the resulting sums, and call the total M(s, oti, . . . , a 2k ). 

(6) The conjecture is 

/oo /*oo 
Z(§ + it, a 2k )g(t) dt= M{\ + it, a u . . . , a 2fc )(l + 0{r^ e ))g{t) dt, 
-oo J — oo 

for all e > 0, where g is a suitable weight function. In other words, L(s, a) and 
M(s, a) have the same expected value if averaged over a sufficiently large range. 

Notes: 

(1) In order to see the structure of these mean values, it is necessary to include the 
shifts ctj. One can obtain the moments of L(| + it) by allowing the shifts to tend 
to 0. Because of the shifts aj we avoid higher-order poles in our expressions. 

(2) The recipe applies to either the L-function or the Z-function, and we give examples 
of both cases. The Z-function case can be directly obtained from the L-function, 
although the reverse is not true in general. 

(3) For the approximate functional equations in the recipe, one can ignore the range of 
summation because it will just be extended to infinity in the final step. 
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(4) We do not define what is meant by a "suitable weight function" , but it is acceptable 
to take g(t) = gr{t) = f(t/T) for a fixed integrable function /. In particular, one 
can take / to be the characteristic function of the interval [0, 1], obtaining the mean 
value J Q T Z(| + it, a)dt. From this one can recover a fairly general weighted integral 
by partial integration. 

(5) The error term 0{t~^ +e ) fits with known examples and numerical evidence. See 
Section [SJ 

(6) The above procedure is a recipe for conjecturing all of the main terms in the mean 
value of an L-function. It is not a heuristic, and the steps cannot be justified. In 
particular, some steps can throw away terms which are the same size as the main 
term, and other steps add main terms back in. Our conjecture is that all of those 
errors cancel. 

2.2. Moments of the Riemann (-function. We illustrate our recipe in the case of the 

Riemann zeta-f unction. In this section we consider the most familiar case of moments of 
C(| + ^0- m Section |2~31 we relate this to moments of Z{\ + it) and repeat the calculation 
for the Z-function of an arbitrary primitive L-function. 

Consider 

(2.2.1) Z(s, a) = ((s + «!)••• ((s + a k )((l - s - a k+1 ) ■ • • ((1 - s - a 2k ), 

where a = (cki, . . . , ac2 k )- Note that this is slightly different notation than given in (|2.1.1|) . 
Our goal is a formula for 

/oo 
Z{\ + it ) a)g{t)dt. 
-oo 

For each (-function we use the approximate functional equation 

(2.2.3) C{s) = ^ ^7 + X(s) ^2 ^7 + remainder. 

m n 

Recall that we ignore the remainder term and the limits on the sums. Multiplying out the 
resulting expression we obtain 2 2k terms, and the recipe tells us to keep those terms in which 
the product of y-factors is not oscillating rapidly. 

If s — z + it with z bounded (but not necessarily real) then 

(2.2.4) x (s) 
and 

(2.2.5) X (l " s 

as t — > +oo. We use the above formulas to determine which products of x( s ) an d xO- ~~ s ) 
are oscillating. 

One term which does not have an oscillating factor is the one where we use the "first 
part" of each approximate functional equation, for it does not have any y-factors. With 



2n 



e lt+ ™'\l + 0{\)) 



2tt 



^- m/4 (l + 0(i)), 
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s = | + it, that term is 

rni s ~ ai ■ ■ ■ m- k s - ak n- l+ak+1 ■ ■ ■ n^ 1+£l2t 



m 1 ,...,m k 
n lt ...,n k 



(2.2.6) = m J 



-i-Qfc -!+a fc+ i -\+<x 2k nx---nk 
' m k n i " ' n k 



m x ,...,m k 
n\,...,n k 



m 1 ---m k 



According to the recipe we keep the diagonal from the above sum, which is 
(2.2.7) 



-4— <*1 -\~ a k -5+Ofe+l -|+"2fe 

//ti ...hi 11- ...11 

wn— m k =m— n k 



If we define 

(2.2.8) R(s;a)= £ 



m\ +ai ■ ■ ■ mt +ak n s - ak+1 ■ ■ ■ n s ~ a2k ' 

where the sum is over all positive m 1; . . . , m^, n 1; . . . , n k such that mi • • • = ni • • - rife, 
then i?(|;a) is the first piece which we have identified as contributing to the mean value. 
(The sum in equation ()2.2.8|) does not converge for s — |. See Theorem 12.4. II for its analytic 
continuation.) 

Note that the variable s in equation 1)2.2.8)1 should not be viewed the same as the variable 
s — n + it from the previous equations. We are employing a trick of beginning with an 
expression involving s and 1 — s, noting that we will later be setting s — |, so instead we 
consider an expression only involving s, which later will be set equal to |. This same trick 
will appear in Section 14.11 when we consider more general mean values. 

Now consider one of the other terms, say the one where we use the second part of the 
approximate functional equation from ((s + a\) and the second part from £(1 — s — a k+1 ). 
By flga and (12~?3D . 

(, ^ ~a 1 +a k+1 
77- 
ZTT 

which is not rapidly oscillating. Using this and proceeding as above, the contribution from 
this term will be 

— j R{\] ttfc+l, 0(2,..., CKjfc, «i, a k+2 , • • • , « 2 fc). 

More generally, note that 
(2.2.11) 

/ . s-i(j-K)t / j. \ - E Pj+Y. ij 

X(s+(3i) ■ ■ ■ x(s+PjMIs-Ji) ■ ■ ■ x(1-«-7jc) ~ ( «- ) e*( J -^/ 4 - 

\z7ie j \2tt 

which is rapidly oscillating (because of the it in the exponent) unless J = K. Thus, 
the recipe tell us to keep those terms which involve an equal number of xi s + a j) an d 
x(l — s — a k+ j) factors. This gives a total of ( 2 fc fc ) = ^ J=0 Q) terms in the final answer. 

We now describe a typical term of the conjectural formula. First note that the function 
R(s; «i, . . . , a k , a k+ i, . . . , a 2k ) is symmetric in a%, . . . , a k and in a k+ i, . . . , a 2k , so we can 
rearrange the entries so that the first k are in increasing order, as are the last k. Thus, 
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the final result will be a sum of terms indexed by the ( 2 fc fc ) permutations a G S 2 k such that 
< • • • < a{k) and a{k + 1) < ■ • • < a{2k). We denote the set of such permutations 
by S. Second, note that the product of an equal number of x( s + a j) an d x(l — s — ock+j), 
as in ()2.2.9j) . can be written as 

a k +a k+1 -\ ha 2 k) / y- \ |( Q cr(l)H ^ a a(k)- a a(k+i) a cr(2k)) 

(2.2.12) 



27T J \27T 

For example, ()2.2.9j) is the case cx(l) = k + 1, a(fc + 1) = 1, and a(j) = j otherwise. 
If we set 

(2.2.13) VK (z, a, a) = ^— J i£(ar; a^j), . . . , a ff ( 2fc )), 
for 2; = x + with x and y real, then combining all terms we have 

(2.2.14) M(z; a) := (J- J ^ W(z, a, a). 
The recipe has produced the conjecture 

/oo roo 
Z(l + it,a)g(t)dt= / M(§ + it, a)(l + 0{r? +e ))g(t) dt, 
00 */ —00 

with Z(s, a) given in ()2.2.1)1 and M(s; a) given above. 

Note that the exponent of (t/2ir) in ()2.2.13|) has half the atj with + sign and the other 
half with — sign, and the same holds for R(s,a). This allows an alternate interpretation of 
H as the set of ways of choosing k elements from {ai, . . . , a 2k }. 

The general case of Conjecture 11.5. H is stated in terms of the Z-function. We can recover 
the mean value of the Z-function directly from that of the L-function (in this case, the (- 
function). By the functional equation and (|2.2.4|) we see that 

. \ |( a iH hatk—atk+i a 2k ) 



Z(s + ai )---Z(s + a 2k ) = (-1 (1 + 0(1/*)) 

(2.2.16) xC(s + «!)•■• C(s + a k )C(l ~ * ~ ajfc+i) ■ ■ ■ C(l - s - a 2k ). 

f + \ hctfe-Qfe+i c*2fc) 

The factor I — ] can be absorbed into the weight function g(t), so 



2tt, 

we obtain the conjecture 



/oo roo 
Z(s + ai)--- Z(s + a 2k )g(t) dt = ^ W(s, a, a)(l + 0(H+ £ ))3(t) dt, 
-00 «/— 00 

where s = | + it. 

In the next subsection we directly obtain the above conjecture for the Z-function of 
a general primitive L-function, and in the remainder of this section we perform various 
manipulations to put these in the form of Conjecture 11.5.11 
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2.3. Moments of a primitive L-function. Consider the primitive L-function 

(2.3.1) £M=E£=IK(=). (<t>1) ' 

n=l p ^ ' 

We assume a functional equation of the special form £l(s) = 7l(s)L(s) = e£ L (l — s), where 

(2.3.2) 7L(a) = 0*n r (5 a + ^)» 

j'=i 

with {yUj} stable under complex conjugation. Note that we have u>j = |, which is conjectured 
to hold for arithmetic L-functions. We also assume 



X" 



n=0 



(2.3.3) = J](l- 7pj x)- 1 , 

where w is the degree of £ and where |7 PJ | = or 1. Again this is conjectured to hold for 
arithmetic L-functions. 

We are going to evaluate the moments of the Z-function Z c (s) = e~^X(s)~^C(s) : where 



X(s) = 1C{1 ~ S) 



We will have to determine when products of X(s) and X(l — s) are not rapidly oscillating. 
By Stirling's formula 

(2.3.5) r(. + ! t)=e-^-(-) e T(-,)^i-.^_-_ + _j + o^_ 
we obtain 

(2.3.6) *( a ) = Q 1 " 2 * (j) f^-j e »(«+W4) (1 + 0(1)), 

as £ — > +oo. Note that the above expression can be simplified because we have assumed 

E %h) = o. 

Now we are ready to produce a conjecture for 

/oo 
Z c {s + «!)••• Z £ (s + a 2k )g(t) dt. 
-oo 

with s = | + i£. 

By the definition of Z, 

(2.3.8) J] Z £ ( S + a,-) = J] + a*)"* II £ ( s + <*>)■ 
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According to the recipe, we replace each C(s + it) by its approximate functional equation 
and multiply out the product obtaining 2 2k terms. A typical term is a product of 2k sums 
arising from either the first piece or the second piece of the approximate functional equation. 
Consider a term where we have i factors from the first piece of the approximate functional 
equation and 2k — £ factors from the second piece. To take one specific example, suppose it 
is the first i factors where we choose the first piece of the approximate functional equation, 
and the last 2k — I factors where we take the second piece: 

e-**(| + a 1+ it)-* ■ ■ ■ X{\ + a t + it)' 1 * £ ■ ■ ■ £ ^ 



2+ae+it 
ri\ Tl-y 7i£ Tip 



e^X(\ + a l+1 + + a 2k + it)* E i-TU " ' ' E ~ 



a n 2 k 



j._ , { X{\ + a 1 + it)--- X{\ + a e + it) \ » 



2 a e+1 ~it Z-^ ^-a 2k ~it 
n e+ i n £+1 n 2k fl 2k 



X{\ + a w + it) ■ ■ ■ X{\ + a 2k + it) 

>2 3 g\ x . . . a ni - - - a>nfin t+x ' ' ' a n k / • • • Uj \ 

^ ^ A +ai . . . . . A~ a 2k \n e+ i ■ ■ -n 2k ) 

ni n 2k Hi tig ">e+l a 2k 

The recipe tells us to retain only the expressions of this sort where the factor involving X 
is not oscillating. By ()2.3.ti|) the requirement is that £ = k (and in particular 2k has to be 
even), and we have 
(2.3.10) 

X(± + ai + it)...X(l + a k +it) \ 2 / Q^V " 



X{\ + a k+1 + it)... X{\ + a 2k + it) 



Now the recipe tells us to keep the diagonal from the remaining sums, which in (|2.3.9j) 
is the terms where n\ • ■ ■ ri£ — ri£ + i ■ ■ ■ n 2k . So in the same way as the ^-function case in the 
previous section we let 



' a n k O-n k+1 " ' ' 0"a k 



n s+ai . . . n s+olk n s ~ ak+1 ■ ■ ■ n~ a2k 

"■i---n k =n k+1 ---n 2k 1 k fe+1 2k 



(2.3.11) R(s,a)= E 

and 

t^2_ \ f («<r(l)H l-a<r(fc)-a<T(fc+i) o CT (2fc)) 
Q w y \ 
~2~ f R(x;a a{1) ,...,a a{2k) ), 

for a G H, the set of permutation of {1, ... , 2k} with cr(l) < • ■ ■ < a(k) and a(k + 1) < ■ • ■ < 
a(2k). Then 

(2.3.13) M(z;a) = J2 W ( 

z,ce,a), 

and we arrive at the conjecture 

/CO 
M(\ + it,a)(l + 0(t-^))g(t) dt, 
-co 

which is the same form as we obtained in ()2.2.17|) . 
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We will now examine the expressions in these conjectures in detail, rewriting them in a 
more explicit form. 

2.4. The arithmetic factor in the conjectures. We retain the notation of the previous 
subsection. In particular, C(s) is a primitive L-function having the properties listed at the 
beginning of Section 12*31 and R(s; a) is given in (|2.3.11|) . 

Theorem 2.4.1. Suppose \atj\ < 5 for j = l,...,2k. Then R(s; ax, • • • , «2fc) converges 
absolutely for a > | + 5 and has a meromorphic continuation to a > | + 5 . Furthermore, 

k 

(2.4.1) R(s; ax, ... , a 2 fc) = ((2s + a { - a k+j ) A k (s; ax, ... , a 2fc ) 
where 

(2.4.2) A k (s; a u . . . , a 2k ) = ]J ( J\ (1 - p ~^+^) j Sp ( s; ax, ... , a 2fc ) 

p v,j=i / 



(2.4.3) a 1; ... , a 2fe ) = jf J[ C p C p (J^j dO. 

Proof. We assumed |7 M -| < 1, so we have the Ramanujan bound a n < d w {n) <C n £ . That 
implies absolute convergence of R(s; a) for a > \ + 6 + e. 

Since the coefficients of R(s; a) are multiplicative, as is the condition n\ - ■ -n k = n k +i ■ ■ ■ ri2k, 
we can write R(s; a) as an Euler product: 



Q"n\ ' ' ' Cln k Cl nk+1 ' ' ' Q'nik 



n 1 ---n k =n k+1 ---n 2 k 

n t e 



n s +ai . _ . n s-a 2fc 



J?(s; Qti, . . . , Qt 2 fc) = E 

l-"2fc 



1 "ft 



pei(s+ai) . . .pe 2fc (s-Q 2 fc) 



1 + I"?' 2 E „2 S+ai -a fc+J + E 
ij=l F j=2 

(2.4.4) = 1 [ C(2s + oa- a k+j ) 



Cpj(ax, ■ ■ ■ , a 2 fc) 

p2js + 



fe 



X 



nl .. l2 1 ^ c' (ax,...,a 2fc ) 

i + (M 2 - 1) E p ^- afe+i + E — + 



k 

(2.4.5) = Y\ C(2s + a i -a jfc+i )A fc (s;a 1 ,...,a 2 fc), 

say. Above c pJ and c' pj are just shorthand for the (complicated) coefficients in the Euler 
product. Estimating them trivially and using the fact that |a p | 2 is 1 on average, (which 
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is conjectured to hold for primitive elements of the Selberg class) we find that Ak(s;a) is 
analytic in a neighborhood of a = |. 

Finally, we have 

k 

(2.4.6) A k {s- an, ... , a 2k ) = J] J] (1 - p- 2s ~ a ^ a ^)B p (s; an, ... , a 2k ) 

where 

dp e l ■ m m CLp^k dp e k+l • • • CLp £ 2k 



B p (s; an, . . . , a 2k ) = ^ 

= l e j="l2j = l e k+j 



pe 1 (s+a 1 ) . . . pe 2 k(s-a 2 k) 



/ D ei(s+ai) . . . v e 2k (s-a 2 k) I \^ ^ ^ K ^ J J 

JV> ei,...,e 2 fe F ^ \ 3=1 3=1 J 

n\ k oo fc 

/ TT y -^—Ae ,0) TT T +J y e(-e w g) 

/ „ej(s+aj) ^ J y 11 „e fc+j (s-a fe+j ) V K +J / 



i=l ej =0 J 3=1 
1 * 



as claimed. □ 

To summarize, the conjecture for the general mean value an, . . . ,a 2k ,g) involves 
the function M(s; a), which can be written as 

(2.4.8) M{s; a) = ^ W(s; a CT(1) , . . . , a^)), 

where we have written W(s; a CT (i), . . . , a<r(2fc)) f° r ^(s; a). And by Theorem 12.4. 1| 

(2.4.9) W(s; an, ... , a 2fc ) = I — ^— ) A fc (s; «i, . . . , a 2fc ) JJ C(2s + a* - a fc +j)- 

One can see the above elements in Conjecture 11.5.11 m particular the form of B p (s; a), 
and Ak(z) in that conjecture equals Ak(\,a) given above. The overall structure is slightly 
different because Conjecture 11.5.11 is expressed as a multiple contour integral, as opposed 
to a sum over permutations. In the next subsection we show how to write the sum over 
permutations in a compact form. In the following subsection we return to the functions A k 
and write them in a more explicit form. 

2.5. Concise form of permutation sums. As we have seen, our methods naturally lead 
to an expression involving a sum over permutations. In this section we describe how to write 
those sums in a compact form involving contour integrals. Similar combinatorial sums arise 
from our matrix ensemble calculations, and we have previously stated our main results and 
conjectures in this compact form. 

Note that in both of these Lemmas, the terms in the sum on the left side have singulari- 
ties. However, examining the right side of the formula makes it clear that those singularities 
all cancel. 
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Lemma 2.5.1. Suppose F(a;b) = F(ai, . . . , a k ; b±, . . . , b k ) is a function of 2k variables, 
which is symmetric with respect to the first k variables and also symmetric with respect to 
the second set of k variables. Suppose also that F is regular near (0, . . . , 0). Suppose further 
that f(s) has a simple pole of residue 1 at s = but is otherwise analytic in a neighborhood 
about s = 0. Let 

k k 

(2.5.1) K( ai , ...,a k -h,... b k ) = F(a u b k ) J[ J[ /(a, - b 3 ). 

i=l j=l 

If for all 1 < i,j < k, ati — a k+ j is contained in the region of analyticity of f(s) then 

/ J K{0£ a {l), CKo-(fc); «cr(fe+l) • • • a<r(2fc)) = 

o-es 

(2.5.2) Lg! 1 / . . . / ^•••,^ afc - ; ^)A(z 1 ,...^) a ^ fc; 

where one integrates about small circles enclosing the ctj 's, and where H i/ie se£ o/ ) 
permutations a G S^fc smc/i t/iat er(l) < ■ • ■ < cr(A;) and a(fc + 1) < • ■ ■ < a{2k). 

The above Lemma applies to the Unitary case, which has been the subject of this section. 
The next Lemma is useful in the Symplectic and Orthogonal cases, which will be addressed 
beginning in Section I4~4l 

Lemma 2.5.2. Suppose F is a symmetric function of k variables, regular near (0, . . . ,0), 
and f(s) has a simple pole of residue 1 at s = and is otherwise analytic in a neighborhood 
of s = 0, and let 

(2.5.3) K(ax,...,a k ) = F(a 1 ,...,a k ) \ [ /(a* + 

i<«<j<fc 

or 

(2.5.4) K(a u . . . ,o fc ) = F(a 1 , . . . ,a k ) ] [ /(oj + flj-)- 

l<i<j'<fc 

//a, + Oj are contained in the region of analyticity of f(s) then 



^ K(exai,...,e k a k ) 



(_l)k(k-l)/2 2 k f r A(zl... 1 Z 2 k ) 2 UU Z i 

(2.5.5) 1 i - K(z u . . . , z k ) — k ''I k> 11 - 1 J - dat • • • dz k 



and 



(JJej)lf(eiai, . . . ,e fc a fc ) 

e,-=±l i=l 



2 * r r A(z 2 , . . . , ^) 2 n * a 



where the path of integration encloses the ±Oj 7 s. 

The proofs of the lemmas come from the following 
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Lemma 2.5.3. Suppose that F(a; b) = F(a±, . . . , a m ; b%, ... , b n ) is symmetric in the a vari- 
ables and in the b variables and is regular near (0, . . . , 0). Suppose f(s) = - + c + • • • and 
let 

m n 

G(ai, ...,a m ',bi,...b n ) = F(a 1 , ........ b n ) J"J \~\ /( a * ~~ h i)- 

Let E m ,n be as defined above. Then 



(— 1) -sr-^ p G(zi, . . . , z m+n )A(zi, . . . , z m+n ) 

m \ n \ 2-^1 ' KeS ( 2: iv,2m+n)=(a (T (i),...,o (T ( m+ „)) i-rm+n TJ m + n ( z . — a .) ' 

cren m+n LLi=l LLj=l V * j) 



Proof. It suffices to prove that 



\2 



Dp, A(zi, ■ ■ ■ , Z m+n ) I — )'"-" 

r *-et>( 2lv .. i2 : m+n ) = ( Q , (T(1)) ... jQ , CT(m+n) ) m+n m+n ^ ^ \ 1 ) 

lli=l lll = l \ Zi a i 



Lj=l 

since each such term will appear m\n\ times. Consider the case where a is the identity 
permutation. Then the residue is 

\2 



-\ \ m+n. 



the answer will be the same for any permutation a. 

The residue above can be expressed as (27ri)~ m ~ n times an m + n fold integral, each path 
of which encircles all of the poles of the integrand; note that the value of such an integral may 
be calculated by summing the residues and note that there is no singularity when Zj = Zk 
because of the factor (zk — Zj) 2 in the numerator. 

□ 

To obtain the form of Conjecture 11.5.11 from the formulas at the end of Sections 12. HI 
and EH apply Lemma f2.5.1l with K(ai, . . . , ai2k) = W(§> &i, ■ ■ ■ > «2fe)- That is, f(z) = 
((1 + z) and 

F(ai, . . . , a 2 fe) = I —j- I A fc (|; ai,..., a 2 fe) 

(2.5.7) = exp ^wlog ^^y-J • ^XJ a J ~ a k+jj A k {\\ a 2k ). 
We arrive at the general case of Conjecture II. 5. l| 

Conjecture 2.5.4. Suppose C(s) is a primitive L-function having the properties listed at the 
beginning of Section \2. 31 and the mean value or, . . . , a^, g) is given in \2.3.1 ). Then 

(2.5.8) a 2fe ,<?) = J P k (w\og(^Ya)(l + 0(t-^ £ ))g(t)dt, 
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where P k (x,a) and G(zi, . . . , Z 2 k) ar & o,s stated in Conjecture \1.5.1\ except that A k is the 
Euler product 



[2.5.9) ^)=nnnf 1 -^-) / 



"P \ Ij... I P \ I. 

p i=lj=l V r / «v j=1 



r 



Note that for the Riemann ^-function, w = 1 and Q = and C p (x) = (1 — 

so Conjecture 11.5.11 is a special case of the above. Also note that wlog (^p^j is the mean 

density of zeros of C{\ + it), or equivalently the log conductor, as expected. 

It remains to express the arithmetic factor A k in a more explicit form, which we do in 
the next section. 

2.6. Explicit versions of the arithmetic factor. The factor Ak(s, a) in the 2fcth moment 
of a primitive L-function can be expressed in a simple form. 

Recall, see Theorem 12.4.11 that A k is the Euler product 

(2.6.1) A k (s; a) = JJ B p (s; a 2k ) f[f[(l- 2s+ ^ Qk+] ) , 
where 

(2.6.2) B p (s; a,, . . . , a 2k ) = £ f[L p (j^j T p (J^j d9. 
Lemma 2.6.1. If C p (x) = (1 — 7 P x) _1 with \j p \ = 1 then 



U(' 

(2.6.3) »,(,;.„...^)=nn(i- ' ) 



k k / , -'--'-I p2s+Q!j-afc+i 



8=1 j=l N ' m=l i^m 

Corollary 2.6.2. If C p (x) = (1 — 7 P x) _1 wift |7 p | = when p\N and \j p \ = 1 otherwise, 

then 

(2.6.4) 



fc J-l \ r)2s+Qj-afc+i / k k 



m«<*,..., °»> = n e n - ; _ * ii n n 1 



"fc + m J. J. J. J. J. J. X p2s+Qi — Qf fe+ j 

pjTV 771=1 i/m ^ p|AT i=l j=l V ^ 

In particular, if C(s) = L(s, x) with \ a Dirichlet character of conductor N , where the 
Riemann (-function is the case N = 1, then 

( 2 - 6 - 5 ) Ms;a u a 2 ) = ]j(l-^^), 

p\N V y J 

2 2 ( i A 

(2.6.6) A 2 (s; a u a 2 , a 3 , a 4 ) = ((is + ai + a 2 - a 3 - a 4 ) _1 JJ JJ JJ ( 1 2s+ai - a2+J ) > 

p\N i=l j=l ^ F 3 ' 
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and 

A 3 (s; at, . . . , ate) 



J^J^l _ p~Y.l a r- a 3+z QgQl _|_ _|_ pQ3) ^p- Q 4 _|_ p~a 5 _|_ p~«6^ p~4s 

_)_ p- E? ai-a 3+i _|_ p Q 2 _|_p«3^ (j)~ ai -\- p~ a2 -\- p~ a3 ) 

+ {p a4 + p° 5 + p a& ) {p~ ai + p~ as + p~ a& ) - 2^jp- 6s 



3 3 



(2.6.7) xnnni-= 

p[JV «=i i=i v 

For A: > 3 it is not possible to express A k as a finite product of (^-functions. 
Proof of Lemma \2.6.1\ Using C p (x) = (1 — 7 P x) _1 and setting 

(2.6.8) a,- = ^ p and Ofc+j = — — — for 7 = 1,..., k, 
we have 

„i k 

B p (s; ai ,...,a 2k ) = / (1 - e^)?,)" 1 (1 - e^g^)" 1 dO 

Jo j=1 

U j= i Qj Jo J=i 

(2.6.9) = -t^-L/^ff^-l/^)- 1 ^-?^)- 1 ^, 



i=i 



where the path of integration is around the unit circle. Since \qj\ < 1, by the residue theorem 
we have a contribution from the poles at qu+i, ■ ■ • , Q2k, giving 



■ — l) k i 

1 Lj=l m=l i=l j^m 

(2.6.io) = ^n^-^-r 1 n^-^+^-r 1 

m=l i=l i^m 

Since 

A fc A; 

(2.6.H) Yl i 1 ~ Qiik+j) Yli 1 - mk+mT 1 = n JJ(i _ ?<?*+j)> 



A: 
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factoring out 

k 

(2.6.12) Jlil-qiQw)- 1 , 

we have 

(2.6.13) B p (s; a)=(f[(l- mk+j )A £ J[ "f^ " ^ ■ 

\i,j=l J m=lijk m 1k+iq k+m 

Since 

(2.6.14) qj q k+i =p- 2s -^ +a ^ and q k+l q k l +m = p a ^~ a ^ , 

we obtain the formula in the Lemma. □ 

Notice that the special case N = 1, i.e. the Riemann ( function, reads in Corollary 12. 6. 21 



. IK* 

(2.6.15) A k (s;a 1 ,...,a 2k ) = l[Y,n — 



p2s+aj-a k+i 



I — rp&k+i— a k + r. 

p m=l i=£m 

Each local factor 

1 



(2.6.16) A Pjk (s;a) = J2 U 



3=1 



1 1 ( ' p2s+aj-a k+i 



\ — p a k + i — a k + m 

m=l ij^m 

is actually a polynomial in p~ 2s , p~ aj and p ak +^ for j = 1, . . . , k. That this is so in p~ 2s and 
p-aj j g reac [iiy apparent from (J2.6.16)) . The fact that it is also a polynomial in p ak +j follows 
from (j2.4.2j) and (|2.4.3j) . from which 

(2.6.17) A Ptk (s; a x ,...,a k , a k+1 , . . . , a 2k ) = A p , k (s; -a k+1 , . . . , -a 2k , -a>i, . . . , -a k ). 

Setting (3i = —a k+ i,...,P k = —a2k, one has, from the above discussion, that A Pjk is a 
polynomial in , i.e. in p ak +^ for j = 1, . . . , k. Finally, use the fact that if an analytic 
function of several variables is of polynomial growth in each variable separately, then it must 
be a polynomial. 

2.7. Recovering the leading order for moments of (. Conjecture 11.5.11 contains, as a 
special case, a conjecture for the leading order term for the moments of the Riemann zeta 
function. In this section we show that the leading order terms derived from Coniecture ll.5.11 
agree with the leading order terms which have previously been conjectured by other methods. 

As described in Section 11.31 it is conjectured that the mean values of the Riemann zeta 
function take the form 

rT 



[2.1.1) [ \al + it)\ 2k dt = TV k (lo g T)+0(T^ 

Jo 



where Pfc(logT) is a polynomial in logT of degree k 2 . Conrey and Ghosh conjectured that 
the coefficient of the log fc T term is of the form g k a k /k 2 \, where a k is given by ()1.3.2J) . 
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Keating and Snaith used random matrix theory to conjecture that g k is given by (jO.Hjl . 
This leading order term g k a k /k 2 \ will be re-derived here, starting with Conjecture 11.5.11 

Conjecture 11.5.11 implies 

(2.7.2) f T \C(i + it)\ 2k dt= f P fe (log£) dt + 0(T^), 

Jo Jo 

where P k is the polynomial of degree k 2 given by 

/ r r k k 

P k( X ) = u\2(n„A\2k f f A k { Zl , . . . , Z 2k )Y[Y[((l + Zi -Z k+ j) 

J J i=lj=l 



k\ 2 (2m) 



[2.7.3) x ^ 2k - 2k e ~ 2 Z3 ~ Zk+j d *~- dz2k - 



1 lj=l z j 

Our goal is to show that the leading order term of P k {x) is {g k a k /k 2 \)x k2 . 

Using the fact that A k is analytic in a neighborhood of (0, ... , 0) and the (^-function has 
a simple pole at 1 with residue 1, after a change of variables we have 

i — 1 j — 1 

x Aa (* 1 '-"'*») e S}-i dzi ... dz 2k 

(-i) k A k (o,...,o) ( x\k 2 _ x r r A 2 ( Zl ,...,z 2k ) 

= kl2 ^ 2k UJ +0{x )] f"'f (nf =1 n-=i^-w)n?=i# 

xe Ej=i zj-zk+j dzi - ■■ dz 2k 
_ A k (0,...,0) x k 2 , x + q, x -i\\ I J A(z 1 ,...,z 2k )A(z 1 ,...,z k )A(z k+ i,...,z 2k ) 



dO. 



k\ 2 2 k2 (2m) 2k K y "J J I\f =l z 2k 
(2.7 '.4) xe E ^=i Zj ~ Zh +j dz!..., dz 2k . 

Now we need only show that A k (0, . . . 0) = a k and the remaining factors give g k /k 2 \. 
From Conjecture 11.5.11 

( „ 5) A ^u(:-^a-f) k (:-^- k 

For a given p, we concentrate on the integral in the above expression, writing it as a contour 
integral around the unit circle: 

(2.7.6) (-pt/l ^'t'-^ fc 

After expanding the two factors in the numerator around z = p~ 5 and calculating the residue 
we are left with the sum 

(2-") (1 - J)-« £ (* T 7 V - i) f . 



1=0 
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Next one can perform a binomial expansion of (1 — and gather like powers of ^ to obtain 
(2.7.8) 

- «- — ^ I J^~, / h -J- n — m — 1\/ h — 1 \ / h A- m — n — l\\ 

p- m . 



(1 _ iN-2fc+iy^ \^(-l)if k + q ~ m ~ 1 \( k ~ l \fk + m-q-l\ 
P V Q )\k + q-m-l)\ k-1 J 

A simple manipulation of the binomial coefficients and replacing q by m — q gives 



E(-D'( 



k + q — m — 1\ f k — 1 \ / k + m — q — 1 
q J \k + q — m — lj \ k — 1 



q=0 

( „ 9) ^^p-ir-fXT 1 



g=0 
k-l\ 



and this final sum over q is in fact just ( ) (see, for example, [RioJ). Thus, 

(fc-i) 2 fc-i /jU lN a 



(2.7.io) A.(o,...,o)=n(i-i) " eC',; 1 ) 

p V ■ P/ m=0 V 7 



and this is indeed equal to a k defined in (jl.3.21) . 

Now we must identify the remaining terms as gk/k 2 \ defined in (|1.3.3|) . as x — > oo. The 
method applied below as far as (|2.7.14|) follows closely that used for a similar purpose in 
|BHj . Expanding the determinants A(zi, . . . ,Zk) — detfzj 1- )j,m=ii we obtain 

lim^ = — - 1 ^ {■■■{<>■ > > 

a k x k2 (k\) 2 2 k2 (2ni) 2k 



x— »oo 



v 5 J \ Q 

(2.7.11) x ^sgn^z^-.-zfr j ^...^^...cfaafc. 

Here Q and R are permutations of {0, 1, . . . , k—1} and S is a permutation of {0, 1, . . . , 2k— 1}. 

Since the integrand is symmetric amongst z 1; . . . , z fc and also amongst Zk+i, . . . , z 2 k, in 
each term of the sum over Q we permute the variables z±, . . . , Zk so that Zj appears with the 
exponent j — 1, for j = 1, . . . , k. In the sum over S the effect is to redefine the permutations, 
and the additional sign involved with this exactly cancels sgn(Q). We do the same with the 
sum over R, and as a result we are left with k\ 2 copies of the sum over the permutation S: 

lim^ = _i - [../e^^ 

x^oo a^x k 2 k [2-Ki) Ik 



x ^>giU .S).: ; : 2 • • • •-/, 

5 

-(2k-S k ) -(2k-S k+l -l) ^-(2fc-5 2fc -i-(*-l)) / 



-(2fe-5 ) -(2fc-S x -l) -(2fc-5 fc _i-(Jfc-l)) 



x ^+i ^+2 az 1 ---dz 2 k- 
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Since 



(2.7.13) 



W) = 2S / (-*)- e "(-*). 



where the path of integration C starts at +00 on the real axis, circles the origin in the 
counterclockwise direction and returns to the starting point, we can rewrite (|2.7.12jl as 



lim 



Pk{x) 



a k x 



k 2 



^ sgn(5) (v(2k - S )T(2k - S, - 1) • • • T(2k - S^i - (A; - 1)) 

x (-l) s *r(2fe - S fc )(-l) 5fc + 1+1 r(2A; - S k+1 - 1) • • • 
x (-1)' 
1 1 



2 fe2 



2 fe2 



,5 2fc _ 1+fc -l r(2A , _ ^ _ (fc _ 1)} 
1 1 



r(2fc) r(2fc-i) 
1 1 



r(2fc-i) r(2fc-2) 



r(i) 

fk-l 

n 

^=0 



1 

r(o) 



-1 

r(fe+i) r(2fc) r(2fc-i) 

1 -1 1 

r(A) r(2fc-i) r(2fc-2) 



1 



r(2-fe) 



-1 
r(i) 



1 

r(o) 



r(fc+i) 

(-i) fc 
r(fc) 

(-i) 3fc - 2 



r(2-fe) 



(k + £)\ 



(2.7.14) 



(!) 
CY) &) 



-CV) (V) 



:-i)*- x (^J 



The above is a 2/c x 2/c determinant, the first k columns of which are identical to the first k 
columns of the matrix 



(2.7.15) 



(!) 
(V) CV) 



( 2fc °-l) \ 



The matrix f)2.7.15|) is lower triangular and so can easily be seen to have determinant equal 
to one. It is also the inverse of 



(2.7.16) 



/ -Q 
-0 G) 

V-CV) (V) 



(2^—1) 

(£:!) / 



It so happens that matrix ()2.7.16|) has its k first columns identical to columns k + 1 through 
2/c of the matrix in expression ([2.7.14)) . Therefore we can multiply expression (|2.7.14|) by the 
determinant of (|2.7.15|) (which is equal to one) and this simplifies the final k columns of the 
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resulting determinant significantly: 

Pk(x) 



lim 



a k x 



k 2 



■l) k 



n 



2*' \^{k + l)\ 



(?) 
fit 1 ) fit 1 ) 



■1)' 



n 



( ) 

\2k-l) 
\2k-l) 



2 (o) 





fit 1 ) 





(I) 



2 *-i ( *-i) 
2 fe (S) 



2 fc-lg) 



L :) 

\k-l) \ ) 




ffc-1 



(-I) fc -Vl) 



1 
1 



Gil) o 
2( fe !j o 



l)/2 



2 fc 2 



n 



2 2fc-i(2*-i) 2 2fc - 2 ( 2 *- 1 ) ••• 2 fc ( 2 fc A l 1 1 ) 

2 2fe-1^2fc-lj 2 2fc-2^2fc-l-| 



{k+£)\ 



'k-\ 



2 fe (S) 

(V) rr 1 ) 

The matrix above can be decomposed as 



(?) 

\ k-1 ) 

fit) 



2ft") 
2fit) 



I (V) fit') 

/ fit) (V) 
fit) (t 2 ) 

I (?) 



fit. 1 ) \ 

I k-l J 



fit) ) 



fit) \ 
fit) 



/ 



ffl 
(t) 



(?) 





(k-l) / V (-fc+l) (-fc+2) 



l) fe (tx) 



1 





(A) \ 

f fc ) 

) 



The first matrix on the right side is zero in the lower right triangle, and the second matrix 
on the right side is upper triangular. Thus we read that the determinant of the matrix on 
the left hand side is (— l) fc ( fc ^ 1 )/ 2 . Therefore, 



(2.7.17) 



lim 



Pk(x) 

k? 



and this is gk/k 2 \ from (jl.3.3|) . as required. 



fc-i 



n 
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A similar method applies to the orthogonal and symplectic cases. 

3. Families of characters and families of L-functions 

We will describe a particular kind of "family" of primitive L-functions based on the idea 
of twisting a single L-function by a family of "characters." In the next section we provide 
a general recipe for conjecturing the critical mean value of products of L-functions averaged 
over a family and we demonstrate the recipe in several examples. 

Note that we use "character" somewhat more generally than is usually covered by that 
term. 



3.1. Families of primitive characters. We describe sets of arithmetic functions which we 
call "families of characters." 

Let T = {/} be a collection of arithmetic functions f(n), and assume that for each / £ T 
the associated L-function Lf(s) = ^2f{n) n~ s is a primitive L-function with functional 
equation L/(s) = £fXf(s)Lf(l — s) and an Euler product of the form 

(3-1.1) L f (s) = £ ^ = J] IK 1 - / WfT 1 - 

n=l p j=l 

The quantity 

(3-1.2) c(f) = \(e f X f )>(i)\ 

is called the log conductor of /. 

Note that if / = %, a primitive Dirichlet character of conductor q, then the log conductor 



is 



(3.1.3) c( X ) 



logg-log7r + f (1) x even 



logg - logTr + T ( i ) x odd. 

If f(n) = n~ lt then the log conductor is c(n~ %t ) = log ^+0(t~ 1 ). Generally the log conductor 
c(f) scales as the log of the "usual" conductor of /. 

In the case that JF is finite, we require that the data Q, Wj, fj,j in the functional equation 
(ll.l.2j) is the same for all / £ T. In particular, the conductor c(f) is the same for all / £ T. 

In the case that T is infinite, we require that the data Q, Wj, /ij in the functional equation 
(ll.l.2j) are monotonic functions of the conductor c(f). Furthermore, we define the counting 
function M(T) = #{/ : c(f) < T} and require that M(log(T)) = L(T A , logT) + 0(T^ +e ) 
for all e > 0, where A > and F(-, •) is a polynomial. 

If G is a function on JF, then we define the expected value of G by 
(3.1.4) (G(/)> = Jim M(T)- 1 V G(f), 



/-x ' 

o(f)<T 



assuming the limit exists. In the case of a continuous family, the sum is an integral. 
We require that if mi, . . . , are integers then the expected value 



(3.1.5) 8 e (mi,...,m k ) = (f(m 1 )...f(m e )f(m e+l )...f(m k )) 
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exists and is multiplicative. That is, if (mim-2 . . . m k , n\n2 ■ ■ ■ n k ) = 1, then 

(3.1.6) 8e(m 1 n 1 ,m 2 n 2 , ■ ■ .m k n k ) = (5<(mi, . . . ,m k )Se(ni, . . .n k ). 

We sometimes refer to 5 as the "orthogonality relation" of the family. 

The practical use of being multiplicative is that a multiple Dirichlet series with 5e coef- 
ficients factors has an Euler product: 

f 3 l 7 \ Mrai,...,m £ ) = T-r 5 e {p e \. . . ,p ee ) 

mi,..., mi 1 * p ei,...,ei 

We will use the above relation in our calculations. 

To summarize, a family of characters T = {/} is a collection of arithmetic functions, 
each of which are the coefficients of a particular kind of L-function. The characters are par- 
tially ordered by conductor c(/), and the expected values Se(mi, . . . ,m k ) are multiplicative 
functions. 

The following are examples of families of characters, two of which are finite and two are 
infinite. The term "finite family" is somewhat misleading, because those families depend on 
a parameter, and the size of the family grows with the parameter. 

1. The family of t-twists. 

(3.1.8) .F t = {f t (n) = n- ft ,0 <t<T}. 

We have 

1 f T .. [ 1 n = m 

(3-1.9) -/ (m/nfdt = { {m/nr _^ . 

T Jo l Tiog(m/n) otherwise, 

leading to the expected values 

1 n = m 



(3.1.10) (ft(n)f t (m)) = (n-V) = ((m/n) rf ) = 



otherwise. 



Therefore the orthogonality relation is 

(3.1.11) S £ (ni, . . . ,n k ) = 8(ni ■■■n t = n £+1 ■ ■■n k ). 

2. The family of primitive Dirichlet characters. For each positive integer q we set 

(3.1.12) F c h(<?) = {f x ( n ) = x( n ) '■ X is a primitive Dirichlet character mod q}. 
We have 

1 ^* —— Jl nEmmodg and (mn, q) = 1 
(3 ' L13) f E X(n)x(m) = i ^ 

X mod q k 

where the sum is over the primitive characters mod q, and q* is the number terms in the 
sum. This leads to the expected values 

1 n = m and (mn, q) = 1 
otherwise. 



(3.1.14) </ x (n)/xM> = (xHxH) 



Since x( m i)x( m 2) = x( m i m 2) we obtain 

(3.1.15) Se(mi, . . . ,m k ) = 5{mi ■ ■ -m e = m £+1 ■ ■ -m k and (mi • • • m k , q) = l) . 
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Note that the condition in the definition of 5g is not ni\ ■ ■ - mg = m^+i ■ ■ ■ m^ mod q. We are 
computing the expected value as a function of q, so one should think of the rrij as fixed and 
q — > oo. The only way to have mi • • • rri£ = • • • m^ mod q for sufficiently large q is to 
have actual equality. 

Note that by our definition, JF ch (g) is not a family, but it is the union of two families 
consisting of the even characters and the odd characters separately. 

3. The family of real primitive Dirichlet characters. 

(3.1.16) jF d = {fdin) = Xd( n ) '■ Xd is a primitive real character mode?, \d\ < X }, 
where d runs over fundamental discriminants. We have expected values 

(3.1.17) (f d {n)f d {m)) = { X d{n) Xd {m)) = ( Xd (nm)) = j^W 1 + ?) 1 n ™ = U . 

I (J otherwise. 

The calculation in the case nm = □ is nontrivial and was first done by JutilaJJ]. (If one 
were summing over all d then the expected value when nm = □ would be (p(nm)/nm). 

In practice one encounters more restricted families, so we let J r d(+) = {fd '■ d > 0} and 
jF d (— ) = {f d : d < 0}, and also 

(3.1.18) F d {a, N, ±) = {f d E F A {±) :d = a mod N}. 

For the family T^{a, N, ±), evaluating the expected value of Xd{n) can be tricky, so we 
provide some useful asymptotics. 

Below we restrict ourselves to < d < X, but the same asymptotics hold if one restricts 
to -X < d < 0. 

Theorem 3.1.1. Let Q = gcd(a,N) not be divisible by the square of an odd prime. Then 

0<d<x / ^6) <°6 p | 2Ar f 

d— a mod N 

Next, assume further that either N is odd or divisible by at least 8 ( this condition is related 
to the fact that Xrf(2) is periodic mod 8), and say n = gD, with (□, N) = 1, and with all 
prime factors of g being prime factors of N . Then 

(3.1.20) yj* Un) ^ Mg)m ^^^^ Ma , N) u^- 

0<d<x ^ I^CI^C V \2NU F 

d—a mod N 

Here /^(a, N) is determined according to Table \3.f~l[ Consequently, for the family J-'dia, N, ±) ; 

(3.1.21) ( X d(n)} 



otherwise. 



Proof. We first outline the proof for f)3.1.19j) . One can count odd fundamental discriminants 
\d\ < X by using the Dirichlet series 

£>nK)^K vl 

d odd 1 1 p odd V y / S V ' V 



INTEGRAL MOMENTS OF L-FUNCTIONS 



41 



p 


a 


h 2 (a,N) 





a G Z 


3/2 


1 


a = mod 2 


1 




a = 1 mod 2 


2 


2 


a = mod 4 


2 




a = 1 mod 4 


4 




a = 2, 3 mod 4 





3 


a = 0,4 mod 8 


2 




a = 1, 5 mod 8 


4 




a = 2, 3, 6, 7 mod 8 





> 4 


a = 1,5,8,9, 12, 13 mod 16 


4 




otherwise 






Table 3.1.1. The function h2(a,N) that appears in Theorem 13.1.11 where 
N = 2 p N Q , with N odd. 



As in the proof of the prime number theorem, the main contribution comes from the pole at 
s = 1, and one has 



(3.1.23) E* 1 -^*" 

d|<X 
d odd 

Next, assume that N is odd and (a,N) = 1. To count odd fundamental discriminants in 
arithmetic progression, d = a mod N, one imitates Dirichlet's theorem for primes in arith- 
metic progression, looking at linear combinations involving Dirichlet characters mod iV of 

X( d ) TT (i , X(p) 



e ^f= n ( i+ „ 

d odd 1 1 p odd v ^ 

If one wishes to further specify d > or d < 0, one can restrict to |d| = 1 mod 4 or 
|d| = 3 mod 4 respectively, with x ranging over Dirichlet character mod AN. The main 
contribution comes from the trivial character whose corresponding Dirichlet series is 

d odd 1 1 p°dd \ ^ / ; , 2JV \ f 

and whose main pole is at s = 1. Therefore, for iV odd and (a, iV) = 1, we have 

(3 ' L26) E* ^^jvJ^n^TT. 

d = a mod N ri 
d odd 

with the same result for —X < d < 0. 

Next, for N odd and (a, N) = Q > 1, one can write, for d = a mod N, d = diQ. Apply 
the above method to di with < d 1 < X/Q, d\ = (a/Q) mod (N/Q), d\ odd, and, because (i 
is squarefree, the extra condition that (di,Q) = 1. Because of this last condition, the Euler 
product that we need to take in ()3.1.25|) is not just over odd p \ (N/Q) but also p \ Q, i.e. it 
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is still 



p odd 



Hence, if (a, N) = Q, 

X 6 yr p 



( 3 - L28 ) E 1 ~ TJanJckT)^ II 



0<d<X 
d=a mod iV 
d odd 



0(4iv/g) q ^ i 2 i p + r 



and the same result for —X < d < 0. 



Finally, we wish to take into account even d. The set of even fundamental discriminants 
consists of —4 and ±8 times the odd fundamental discriminants. 

Again, assume iV is odd. One can count discriminants, d = a mod N, lying in the interval 
(0, X) by counting odd discriminants lying in (0, X), together with odd discriminants in 
(— X/4,0), (0, X/8) and (— X/8,0). Overall, this gives the same asymptotics as before, but 
with an extra factor of (1 + 1/4 + 2/8) = 3/2. This accounts for line 1 in Table I3~l.ll The 
other lines in the table can be obtained by similar considerations. 

We now apply to obtain (i;il.2()|) and (|3.1.21|) . Consider 

(3.1.29) 



0<d<X 
d~a mod JV 



E 



d=a mod N 



(the following analysis also holds for —X < d < 0). 

Write N = Nl 1 -...- N^, the prime factorization of N, and let g = N^ 1 • . . . • N^ m . Then 



(3.1.30) X d(n) = Xdi^r ■ ■ ■ Xd(N m ) u - X d(D)- 

Now, if Ni is an odd prime, Xd{^i) — Xa(Ni), since d = a mod N, and so d = a mod N{. If 
Ni = 2 we need to be careful because Xd(2) is periodic mod 8. Now we are assuming that if 
iV is even it be a least divisible by 8, i.e. that d = a mod 8, and thus that Xd(2) = Xa(2). 

Therefore 

/o i o-n f\ ,\ rn\ fxa(g) if (d, □) = 1 

(3.1.31) Xd{n) = Xa{g)Xd{a) = < . 

I U otherwise, 

and one gets that (|3.1.29|) equals 

* 

0<d<X 
d— a mod N 

(3.1.32) xM- " -' 



E 



0<d<X 
d=a mod N 
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Since (□, N) = 1, the sum in the numerator can be split into sums d mod NO. Naively, 
one expects to have </>(□) sums, one for each residue class (d, □) = 1. However, if □ is 
even, then only half of these residue classes, namely those that have d = 1 mod 4, contain 
fundamental discriminants, so one only gets 0(D) /2 sums. We thus consider the case that 
□ is odd separately from the case that it is even. Both cases end up giving the same answer. 

Assume that □ is odd. To apply our formula (J3.1.19)) to each of the </>(□) residue classes 
modiVTJ, one needs to compute the various components that go into the formula. 

Given d = a mod N and d = b mod □, one has via the Chinese remainder theorem 
d = a mod NO. Now, Q = (a, N) = (d, N), and (d, □) = 1, so (a, NO) = (d, NO) = Q. 

One also needs to evaluate h 2 (d, NO). Let iV = 2 /3 iVo, with N odd. Now □ is odd, and 
so h 2 (d, NO) only depends on a mod 2^, but this is determined by a mod N. So h 2 (d, NO) = 
h 2 (a,N). Therefore, the numerator of (jrS.1.32j) is asymptotically 

Canceling factors appearing in the asymptotics ()3.1.19j) of the denominator of (|3.1.32j) we 

get 



(3.1.34) Xa(g)H 



P 



p+l 

If □ is even, write □ = 2 A D , with A > 2, and O odd. Now, (d, O) = 1, so d is odd. 
In all cases, according to Table l3~l.ll h 2 (d,NO) is therefore 4. Furthermore, as in the odd 
case, (a, NO) = Q. 

Hence, one gets asymptotically for the numerator of ()3.1.32|) 

(3 - L35) x ' (g) — wmQ)Q^u a j^- 

Since □ is even, N is odd. Hence h 2 (a, N) = 3/2, and the denominator of ()3. 1.32)1 is 
asymptotically 

(3 - L36) ^T^g^n^TT- 

Canceling numerator and denominator, taking special care for powers of 2 appearing in □, 
we get 

(3.1.37) Xa{g) ^l[^. = Xaig) Y[P 

□ 

4. The family of coefficients of holomorphic newforms. 

(3.1.38) F mod (k,q) = {/(n) = X f (n) : ^ n^ 2 X f (n) E H k {q)}, 
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where Hk(q) the set of newforms in Sk(To(q)). A good reference for these functions is 
Iwaniec flwlj . In this family the parameter tending to infinity can be either k, or q, or some 
combination. 

The orthogonality relation here is somewhat subtle, and in fact there are two natural 
ways to average over these characters. In both cases the starting point is the Hecke relation 

(3.1.39) \ f (m)\ f (n) = ^ \ f (mn/d 2 ), 

d\m,d\n 
(d,q)=l 

which imply that any product 

(3.1.40) A / (m 1 )---A / (m fe ) 
can be expressed clS 8b linear combination 

(3.1.41) $>A/(J) 

for some integers hj, and in fact only for j a prime power. Thus, we need only determine the 
expected value of Aj(p'). 

If one averages over Hk(q) in the most straightforward way, then for p \ q, 



(3.1.42) (A,(p*)> 
and more generally, if (n, q) = 1, 

(3.1.43) (A/(n)> = 



p • ? / 2 j even 
j odd, 



n 1/2 n — □ 

otherwise, 



This follows from the Selberg trace formula. However, if one averages with respect to a 
weighting by the Petersson norm: 



(3-1.44) E h * = E */</./>. 



then 

(3.1.45) (A,(p>)> 
and more generally, if (n, g) = 1, 

(3.1.46) (A/(n)> 

This follows from the Petersson formula (see |Iwlj ). if (mn,q) = 1, 

f o i x / \\ / \ u \ i o -fc 5(m,n;cg)J A; _ 1 (47r v ^:/cg) 

(3.1.47) > \f(m)\f(n) = d(m, n) + 27rz > . 

feH k (q) c=l C? 



i j = o 

otherwise, 



1 n = 1 

otherwise. 
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Here Jk is the Bessel function and 

f ma + nc\ 

(3.1.48) S(m,n;c)= e ( J 

ad=l mod c 

is the Kloosterman sum. Since the Petersson weighting leads to a somewhat simpler ex- 
pression, we will consider that weighting in our example. When passing from the Petersson 
formula to the expected value, using the Weil bound for the Kloosterman sum and the fact 
that Jk has a kth order zero at 0, we see that for fixed m, n the sum on the right side of 
f!3.1.47|) vanishes as k — > oo or q — > oo. 

Let 

(3.1.49) S(mx, . . .,m k ) = (A/(mi) • ■■X f (m k )). 

So in the Petersson weighting, 6 (mi, . . . , m k ) is the the coefficient bi of A/(l) = 1 in ()3.1.41j) . 
One can use the Hecke relations to show by induction that 5 is multiplicative in the sense 
of 13.1.61 Thus, we only need to know S on prime powers. 

Lemma 3.1.2. With respect to the Petersson weighting, if p \ q then 

s(p m \..., P ^) = - [\^eT] sin ^ + Q 1)e de 

7r /n - LJ - sine/ 

OP 71 k p i{mj+l)6 _ -i{rrij+l)e 

(3.1.50) = -jf sin^n 6 eU _% d6. 

For the unweighted sum we have 



(3.1.51) S(p^, . . . ,p^) = i r ! TT t 

If p\q then 5(p mi , . . . ,p mk ) = unless mi = ■ ■ ■ = m k = 0. 
Proof. We only give the details for ()3.1.5()j) . Beginning from 

oo 

(3.1.52) C p (x) = J2 X f (p>)x 3 = (1 - e^fx)' 1 (l - e^^x)' 1 

j=0 

we have A/(p J ) = —^g^ 2 = Uj(cos6f, p ) where C/j is the usual Tchebychev polynomial. 
Then 5 (p mi , . . . , p mk ) = Co where 

(3.1.53) U mi U m , 2 . . . U mk = ^ c e U e . 

e>0 

If we evaluate both sides of this equation at cos 9 and integrate from to 7r with respect 
to the measure - sin 2 6 d6, then the result follows from the orthogonality of the Tchebychev 
polynomials with respect to this measure. □ 

5. The family of coefficients of Maass newforms. 

(3.1.54) J~m(q) = {fin) = X f (n) : X f (n)K tR (2n\n\y)e 2mnx G H(q)}, 
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where H(q) the set of Maass newforms on To(q). A good reference for these functions is 
Iwaniec |lw2j . The orthogonality relation is derived from the Kuznetsov trace formula. See 
Chapter 9 of jlw2]. 



3.2. Families of L-functions. We use a family of characters to create a family of L- 
functions in the following manner. 

Begin with a fixed primitive L-function 



OO / . 

a n T-r - / 1 



(3-2.1) £ W = ES = n^(^). (»>!)■ 

We assume that 



n=l p v ^ 



(3.2.2) C p (x) = V s " = II^ 1 _ ^y 1 ' 

n=0 j=l 

where w is the degree of C and where |7 PJ | = or 1. Assume £(s) satisfies functional 
equation 

(3.2.3) C(s)=eX(s)C(l-s), 
as described in Section IT~TI 

We create a family of L-functions by twisting £ by a family of characters. Let T = {/} 
be a family of characters, with the properties described in Section 13.11 The twist of C by / 
is denoted by C(s, f) and is given by a Rankin- Selberg convolution: 

v w oo I „\ 

(3.2.4) a*, f) = n n IK 1 - ^/pt 1 = e ^f 2 - 

p i=l j=l n=l 

Note that if w — 1 or v — 1, as will be the case in our detailed examples, 

fln/(n) 



(3.2.5) ^(a,/) = E 



n=l 



We require that C(s, f) is an L-function. That is, our family of L-functions must consist 
of L-functions! In particular, C(s, f) satisfies a functional equation 

(3.2.6) C(s,f) = e f X f {s)Z(l-s,f), 

as described in Section fl~Tl As part of our definition of "family", we make a restrictive, but 
natural, assumption on Xf. We have Xf(s) = 7/(1 — s)/jf(s) where 

w 

(3.2.7) 7/(s) = g^JJr(i S + /U , i/ ). 

We assume that w is constant, and each of Qj, $t(fi,jj), and jj) is a monotonic function 
of the conductor c(f). In practice this will mean that each of those quantities will either be 
constant or will be tending to infinity with the conductor. 

For example, the collection of all real primitive Dirichlet L-functions L(s, Xd) is not a 
family, because for d > we have nx t d = and for d < we have = 1, so ni t d is not a 
monotonic function of c(xd)- So we consider these two families separately. 
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Finally, we will make use of an approximate functional equation of shape 



(3.2.8) C(s, f)=J2^P- + £ f X f^ E ^TT + remainder. 

Note: we are not claiming that the "remainder" in the above equation is small, nevertheless 
we will ignore the remainder in our calculations. 

4. A RECIPE FOR CONJECTURING MOMENTS, WITH EXAMPLES 

We give a general recipe for conjecturing the moments of a primitive family of L- 
functions, and then apply the recipe to several interesting examples. 

4.1. The general recipe. Suppose £ is an L-function and / is a character with conductor 
c(/), as described in Sectional So 

(4.1.1) Z c (s,f) = e-*X f (s)-*£(s,f), 
which satisfies the functional equation 

(4.1.2) Z £ (s,f) = Z2{l-s,f), 

so Zc(s, f) is real on the 1-line. Note that the square root of e$ involves a choice which must 
be made consistently. 

We consider the moment 

(4.1.3) + ai > f)---Zc{\ + a k , f)g{c{f)) 

where g is a suitable test function. The recipe below also applies to averages of products of 
C{\ + a, /). The sum is an integral when T = Tt- 

Here is a recipe for conjecturing a formula for the above moment: 

(1) Start with a product of k shifted L-functions: 

(4.1.4) Z f (s, ai, ■ ■ ■ , a k ) = Z c (s + a x , f) ■ ■ ■ Z c (s + a k , f). 

As we will demonstrate in our examples, the recipe applies to the Z-function as well 
as the L-function. 

(2) Replace each L-function with the two terms from its approximate functional equation 
(|3.2.8|) . ignoring the remainder term. Multiply out the resulting expression to obtain 
2 h terms. Write those terms as 

(4.1.5) (product of Ef factors) (product of Xj factors) (summand). 

m,...,n k 

(3) Replace each product of ^/-factors by its expected value when averaged over the 
family. 

(4) Replace each summand by its expected value when averaged over the family. 

(5) Complete the resulting sums, and call the total M(s, aci, . . . , a2k)- 

(6) The conjecture is 

(4.1.6) Z fd «i. ■ ■ ■ » <*»)9(c(f)) = M fd «!>•••> «2 fc )(l + 0(e^ + ^))g(c(f)), 
for all e > 0, where g is a suitable weight function. 
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In other words, Zf(s, a) and M/(s, a) have the same value distribution if averaged over 
a sufficiently large portion of the family. Note that the dependence of Mf on / only occurs 
in the product of Xf factors. 

As we mentioned earlier, some of the individual steps in this recipe cannot be rigor- 
ously justified. Only by using the entire recipe does one arrive at a reasonable conjecture. 
In particular, we ignore off-diagonal terms which actually make a contribution. However, 
comparison with examples in the literature, random matrix moments, and numerical data, 
suggests that the various errors in our recipe all cancel. The underlying cause for this remains 
a mystery. 

We will apply the recipe to several examples, but first we do the initial steps of the recipe 
in some generality. 

For each Z^ substitute the expression in (|4.1.1|) . After replacing each £(s,f) by its 
approximate functional equation (jH.2.8)) . multiply out the product. A typical term is a 
product of k sums arising from either the first piece or the second piece of the approximate 
functional equation. Consider a term where we have I factors from the first piece of an 
approximate functional equation and k — i factors from the second piece. To take one 
specific example, suppose it is the first I factors from which we choose the first piece of the 
approximate functional equation, and the last k — £ factors from which we take the second 
piece of the approximate functional equation: 

S(/) sr^ a nAf) 



n\ 1 rig I 



(4.1.7) xej Xf(s + a e+1 )»--'Xf(s + a k )' ^ i-I-a m ' ' ' 2^ ~ 

n e+1 U l+\ n k H k 

Rearranging this expression and using the fact that Xf(s) = Xf(l — s) _1 , we have 



ef'UXfis + aj)^ J] X f (l-s 

3=1 j=l+l 



CXj) * 



/ 41g x x v ani ^ " ' a ^AD^±i Him ank ill 

ni,...,n fe n l ' " n e ,L i+l ' ' ' n k 

A little trick: since we will eventually set s = |, we replace the above expression by 

l k 



ef £ Hx f (s + a^ II X f( s 

j=l i=l+l 



^ I x a m(/) ' "S(/)S+i(/) ' " a nt(/) 



S + Ql S+OLl^S— Q£ + l „S— Ctk 

n u ..,n k U l ■■■ n l U i+1 



It is expression (j4.1.9|) . and the corresponding pieces from the other terms when multiplying 
out the approximate functional equation, which will appear in the final conjecture, evaluated 



at s — 2 • 



INTEGRAL MOMENTS OF L-FUNCTIONS 



4!) 



-—I 

Now consider the product of Ef factors Ef , which according to the recipe should be 

replaced by its expected value. An important issue is the choice of the square root. We 

i 

believe that there is a natural choice of e? so that the following hold. 

— — £ 

• Unitary case: the Ef are uniformly distributed on the unit circle, and (ej ) = 

unless | — i = 0. In particular, k must be even. There will be ( fc , 2 ) terms in the final 
answer. 

• Orthoconal case: Ef is constant (1 or —1) over the family, or Sf = 1 for approximately 

half the / and Ef = — 1 for the other half. We have (eJ ) = unless | — £ is even. 

In particular, k must be even and there will be 2 k ~ 1 terms in the final answer. 

— — £ 

• Symplectic case: Ef = 1 for all /, and (e? ) = 1 for all k and L There is no 
restriction and there will be 2 k terms in the final answer. 

Note that if we are considering the L-function, instead of the ^-function, then the issue 

i 

of Ef does not arise and the calculation is somewhat easier. See ()4.5.4|) and the discussion 
following. Also note that in the Unitary and Orthogonal cases, odd powers of the Z-function 
will average to zero, while odd powers of the L-function will not. 

The recipe now says to replace the summand by its expected value when averaged over 
the family. That is, we replace 



(4-1.10) a ni (f) . ..a ne (f)a ne+1 (f) ■ • -o» fc (/) 

by its expected value when averaged over the family. In practice, this will be of the form 

(4.1.11) c(^)5 £ (n l7 ...,n k ) 

where c(T) depends only on the family, and where the 8e are multiplicative functions, i.e. 

(4.1.12) 8 e (mini, m k n k ) = 8 E {m u m k )8 t {n u . . . , n h ) 

whenever (mi • • • m k , n\ • • • n k ) = 1. 

The final step is to extend the range of summation. This produces one term in the 
conjecture. By considering the other terms when multiplying out the approximate functional 
equations, one arrives at a conjecture for the original mean value. 

Although the above steps have produced an answer, it is not written in a particularly 
usable form. There are three more steps to put the conjecture in the form of Coniecture ll.5.11 
writing the main terms as an Euler product, identifying the polar part, and expressing the 
combinatorial sum as a multiple integral. 

Since the 8e are multiplicative, we can write the main term as an Euler product. Specif- 
ically, 

/ 4 1 13) >p 8 e (n 1 ,...,n k ) = -j-r 8 e (p e \ . . . ,p £fc ) 

V ' ' ' 2-^ n S+ai ■ ■ ■ n S+ ° k ^ pei(s+a 1 )+-+e k (s+a k ) 

m,...,n k =l 1 k p ei,...,e k =0 F 

assuming that Jfe is sufficiently large. 

Next we determine the leading order poles. It usually turns out that 8i(p £l , . . . ,p £k ) = 
when Yl e j = 1- Thus, the first poles come from those terms where Y2 e j = 2- This happens 
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in two ways, so the rightmost poles are the same as the poles of 

(«■") n n('^)*nn('^). 



l<i<j<k p v ' j=l 

In practice the first factor has simple poles at ~ — + acj), and the second factor is either 

.,. ./. v..- l i ,„», a simple pc^ 



regular in a neighborhood of a = ~, or else it has a simple pole at s = | — a,-. Accordingly, 



we factor out either 
(4.1.15) JJ C(2s + a< + «i) or ] ) ((2s + a { + ay). 

l<i<j<k l<i<3<k 

The remainder is the in our conjectures, and it is regular in a neighborhood of a = |. 

Having identified the polar part of our main terms, we can apply the lemmas in Sec- 
tion E31 to express the sum of terms as a contour integral. The result is an expression similar 
to Conjecture 11.5. II 

We have already seen this procedure in Section 12.21 for the case of mean values of the 
zeta-function. In the following sections we carry out example calculations for families of each 
of the three symmetry types. 

4.2. Unitary: moments of a primitive L-function. The recipe for mean values in Sec- 
tion 12.11 is a special case of the general recipe. To see this, note that if f t G Tt then 
ft{n) = n~ lt , so C(s, ft) = C(s + it). From the functional equation C(s) = eX(s)C(l — s) we 
obtain the functional equation 

(4.2.1) C{s,f t )=e t X t (s)Z(l-s,f t ), 
where 

(4.2.2) e t = eX(l + zt), and X t (s) = ||±f| - 

Note that these satisfy the requirements \X t {\ + iy)\ — 1 for y real, with X t {\) = 1 and \e t \ = 
1. Also note that the log conductor of C(s, f t ), defined as | (e i: A' t ) / (^) | , equals + it)\, in 
agreement with the usual notion of conductor in t-aspect. 

Replacing the product of e r factors by their expected value is the same "keep the terms 
where the product of the y-factors is not oscillating." Thus, after multiplying out the 
approximate functional equations there will be ( 2 fc fc ) terms which contribute. In each of those 
terms replacing the summand by its expected value is the same as "keeping the diagonal." 
Thus, we arrive at the same conjecture as before. 

4.3. Unitary: all Dirichlet L-functions. We apply our recipe to conjecture the average 



(4.3.1) Z x(b®u 



X mod q 
X even or odd 



where the sum is over either the even or the odd primitive Dirichlet characters mod q and 

(4.3.2) Z x (s\ ax,..., a 2 k) = Z(s + «i, x) • • • Z{s + a 2 k, x)- 

Here Z(s + ctk,x) — ( £ x^x( s )) ~^L{s, x) where L(s,x) = £ xX x ( s )L(l — s,x)- Note that 
£ x = r (x)<?~ 5 ? which is uniformly distributed on the unit circle. 
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Following the general discussion in Section UHJ equation ()4.1.9|) specializes in this case to 

l 2k 

(4.3.3) e k - 1 J] X x {\ + a,-)"* J] X x (§ - a,)"* £ j xK)-xK). 

3=1 j=£+l ni,...,n 2h n * . . . 

According to the recipe, we replace e£ - ^ by its expected value. Since the e x are uniformly 
distributed on the unit circle, the expected value is 1 if £ — k and otherwise, so we keep 
( 2 fc fc ) terms. 

Next we replace the summand by its expected value, which is 



(4.3.4) 



5(n u . . . ,n 2k ) = (x(ni)---x(n jk )x(n fc+ i)---x(n 2 fe)) 

1 n 1 ---n k =n k+ i---n 2 k and (ni • • • n 2 k, q) = 1 
otherwise. 



The above is almost identical to the conjectures obtained for the mean values, in t-aspect, 
for a primitive L-function. So one obtains the same formulas as appear in Conjectures 11.5. II 
and 12.5.41 the only changes being that one omits the factors p\q in the Euler product A k) 
and one must use the factors X x (^ ± <x,-)~2. Specifically, in Conjectures 11.5.11 and l2~5~4l 
a simplification occurred by use of equations ()2.2.4|) and ()2.3.6|) . If those conjectures are 
written in terms of n^(| ^ z j)~ 1 ■> then the Dirichlet L-function moment conjecture would 
be obtained by substituting with F] X x {\ ±Zj)~^ Note that we are considering the averages 
over the even and odd primitive characters separately, so in the sum X x only depends on the 
conductor of x- See the comments following the Theorems in Section [T~o1 for more discussion 
on these X-factors and conductors. 

4.4. Symplectic and Orthogonal: quadratic twists of a real L-function. Next we 
consider what happens when we average the shifts of central values of C(s) twisted by the 



family of quadratic characters XdM = I — ] , d < a fundamental discriminant. Here 

w 

Xd(n) = ^— j is the Kronecker symbol which is a primitive Dirichlet character of conductor 

We will see that the family can be either Symplectic or Orthogonal, depending on the 
particular L-function we start with. 

Again, let 

oo 

(4.4.i) ^)=E!i=n^( i /p s ) 

be a primitive L-function and note that 



n=l 



(4.4.2) C(s, Xd ) = ^ n{ ~ n) 



oo 

n/ 



n=l 



We assume that £ is real, i.e. £ = C, as this case is relatively easy to deal with from a 
fairly general perspective. Thus, 

(4.4.3) C(s) = eX{s)C{l - s), 
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where e = ±1. The twisted L- function is expected to satisfy a functional equation of the 
form 

(4.4.4) C(s, X d) = e d X d (s)C(l - s, X d). 
It is further expected that 

(4.4.5) X d (s) = \d\ w ^X(s,d), 

where there are only finitely many possibilities for X(s,d). By our definition of "family" 
we require that the parameters in the functional equation be monotonic functions of the 
conductor. Since there are only finitely many choices for X(s, d), we must restrict to averages 
over sets of d for which X(s,d) is constant (as a function of d). In the situation described 
here, it is believed that there exists an integer N, depending on C, such that e d and X(s, d) 
only depend on the sign of d and on (d mod N). Thus, we will consider the averages 

(4.4.6) L d faa 1 ,...,a k )g(\d\), 



d<0 
d = a mod N 



(the following analysis holds also for d > 0) where denotes a sum over fundamental 
discriminants d, and 

(4.4.7) L d (s; ai, . . . , a k ) = Z c (s + a u Xd) ■ ■ ■ Z c {s + a k , Xd)- 

Note that e d = e a , which may depend on the sign of d. If N is even we are insisting further 
that it be divisible by at least 8. 

Following the general discussion in Section |4~T1 equation (|4. 1.9(1 specializes in this case to 
t k 

(4.4.8) ef £ l[Ms + a^ ]J X d (s - a 3 )^ J^'.'.'^ XdM ■ • - X d(n fc ). 

j=l j=£+l ni,...,rife 1 k 

-—£ 

According to the recipe, we replace ej by its expected value. We have assumed (by our 

choice of a mod N) that e d = e a for all d, so the expected value is el and we will have 
a contribution from all 2 k terms. (That expression is more transparent is you separately 
consider the cases e a — 1 and e a = —1). 

The next step in the recipe is to replace the summand by its expected value. Since 
Xd( n i) ' ' " Xd( n fc) — Xd{n>i ■ ■ • n k ), from equation 1.21(1 we have the expected value 

Xa(g) n p p( 1 + n l ---n k = gU 

otherwise, 



(4.4.9) (Xd{ni)---Xd{n k )) 



where (N, □) = 1, and all the prime factors of g also being prime factors of N. So the 
contribution from the term where we use the first part of the approximate functional equation 
for the first I factors, and the second part for the rest, is 

e k 

k—£ -| — r 1 -| — r 1 

(4.4.10) el YY X d {s + aj)~* | j X d (s - 5 i2 fciV ( S ; ai , . . . , a e , -at+i, ■ ■ ■ , -ac k ) 

j=l 3=1+1 

where 

(4.4.11) R k , N (s;a u ...,a k ) = J2xa( 9 ) ^ JlC.^ IK 1 + 
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Adding up all 2 k terms we obtain 

k 

(4.4.12) M(s; a 1 ,...,a k )= ^ si 9 n ({ e j}) JJ x d{\ + Cj a j)~* ^fe,Jv(s; • • • , e k a k ), 

€i=±l j = i 

where 

' 1 e„ = 1 



(_1)*E* £o = -l. 



(4.4.13) sign{{ej}) ■■ 
So the recipe has produced the conjecture 

(4.4.14) L d (±,a 1 ,...,a k )g(d)= M(±; ai , . . . ,a k )(l + 0(\d\-^))g(d). 



d<0 d<0 
d=a mod N d=a mod N 



To put the conjecture in a more useful form, we now write R k ,N as an Euler product, 
and then express the main term as a contour integral. 

We have R k ,N = Yl p Rk,N,p, which naturally separates into a product over the primes 
which divide N and a product over the primes which do not divide N . The p-factor when 
p\ N is 

w»> = (i+(i+i) _1 f: e n^&r 

V V F/ 3=1 ei+...e fc =2j i=l F 



1 + 



' ^ 1 J~J 



^/ \^ j=0 ei+...e k =2j i=l ^ 



) - H)~'(H(n^) + n^))) 

Similarly, the p-factor when p|iV is 

(4.4.16) ^ = n^(^)- 



The above expression will enable us to locate the leading poles of R k> N- Consider the 
expansion of R k ,N, P (for p \ N) in powers of 1/p. The expansion is of the form 



i + E^+ E 3^k+o(p-'-^)+o ( p-^) 

j=l y l<i<j<k 

( 4 ' 4 ' i7 > =n( i +^)x n ( i +55^)><(i+o( P - i - 2 " + ')+o(p--«)) 

7 = 1 \ r / l<i<j< k \ r / 



We assume that 

(4.4.18) nf 1 



[a pJ 



p" 
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has a simple pole at s = 1. This is conjectured to be equivalent to C(s) being a primitive L- 
function, and this is the key place where the assumption of primitivity enters the calculation. 
We also assume that 



(4.4.19) U( l + a f 

P v 

has a pole of order <5 = 0orlats=l. 

In general, 5 is expected to be or 1 according to whether the symmetric square L- 
function of C(s) is analytic at s = 1 or has a simple pole at s = 1. If C(s) is a degree 1 
L-function, (that is, the Riemann (^-function or a Dirichlet L-function), then 5 = 1. If C(s) 
is associated to a GL(2) automorphic form, then 5 = in general (except possibly when L 
is a dihedral Artin L-function associated to an weight 1 modular form). 

Note that W v (l + 0{p~ 1 ^ 2s ) + 0(p~ 3s )) is regular for a > |. Thus the total order of the 
pole of the above product at s = | when ati = ■ ■ ■ = = is \k(k — 1) + 5k. Accordingly, 
we factor out appropriate zeta-factors and write the above product as 

RkMs) = n cps+oi+aAiiRwM n f 1 - V 2 S+ i +aj ) 

l<i<j<k p l<i<j<k ^ ' 

if 5 = 0, and as 

RkAs) = II C(2s + a i + a j )l[R k , N , p (s) II f 1 ~ ) 

if 5 = 1. In the first case above the family is Orthogonal, and in the second case it is 
Symplectic. 

In summary, we are led to conjecture 
^2 z d\ + on, Xd) ■ ■ ■ Z c {\ + oi k , Xd)g(\d\) 

d<0 
d=a mod N 

k # 
= ^2 sign({ei}) Y[X(± + ejU^a)'* ^ R k,N{\, ei«i, . . . , e fc afc)|d|^= ie ^ 

ej =±l i=l d<o 

d— a mod iV 

x(l + 0(|d|-^))»(|d|). 

(4.4.22) 

The analogous sum over d > leads to a similar conjecture. Here si<7n({ej}) is given in 
(14.4. 13(1 and in either case we can use Lemma f2. 5. 21 to write the sum as a contour integral. 

In the case that C(s) is the Riemann zeta-f unction, the above reduces to Conjecture ll.5.31 



4.5. Orthogonal: L-functions associated with cusp forms. Recall that the set of prim- 
itive newforms / G S n (T (q)) is denoted by H n (q). In this section we consider the shifted 
moments of Lf(s) = Y^=i A/(n)n~ s near the critical point s = | averaged over / G H n (q). 
Note that in the language of Section l3~2l these L-functions are the twists of £(s) by the family 
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of characters H k (q) and would be denoted as ((s,f). However, we use the more common 
notation here. 

The functional equation is 

(4.5.1) L f (s)=e f X(s)L f (l-s), 

where £/ = — y/q\f(q) = ±1. 

We consider the "harmonic average" 

(4.5.2) Y, h + + 

which attaches a weight (/, to each summand. That is, 
(4-5.3) E h * = E */</'/>■ 

feH n (q) feH n (q) 

Following the general discussion in Section I4.lt equation (|4.1.9|) specializes in this case to 



k 



(4.5.4) e)-" [] X(s - a,)' 1 £ 



A/(ni) • ■■\f{n k ) 



i=M-l n lv ..,n fe n l "' n fc 



According to the recipe, we replace e'j 1 by its expected value. Since Sf is randomly ±1, the 
expected value is unless k — £ is even. Thus, we will have 2 fc_1 terms in the final answer. 

Next we replace Aj(ra 1 ) ■ • • Xf(n k ) by its expected value. This is given in Lemma 13.1.21 
After factoring into an Euler product and summing the relevant geometric series we see that 
()4.5.4|) contributes Yij=i X(s — ctj)~2 R(a\, . . . , ag, — c^+i, . . . , — a k ), where 

R(s,ax, ...,a k ) 

k w k e i8 ( 1 fgj) ~* - e -i8 (l_ e-' e \ ~ 1 

(4.5.5) =n^ + %)- 4 n^/ ^ e-° d " 



Here remember that s will eventually be set to | and A(s) = A(l — s) x . Adding up all 
2«-i ^ erms we obtain 

(4.5.6) M(s;ai,...,a fc ) = JJX(s-a.,-)~2 ^ _R( S; eiQl) . . . ; e k a k ), 

j=l S=±! 

so the recipe has produced the conjecture 

J] h L/d + aO.-.L/d + afc) = M{\ ] a ll ... ) a k ){l + 0{nq)-^ +£ ) 

feH n {g) f€H n (q) 

(4.5.7) = (l + 0(ng)^ +£ ) M(|; ax,..., a fc ) 
Summarizing, we have 
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Conjecture 4.5.1. With A k (a±, . . . , a k ) as in Conjecture M.5. we have 
E h L f {\ + a l )...L f (\ + a k ) = 

feH n (q) 

k k 

j=i Cj =±i j= i 

(4.5.8) x ] ] ({l + E i ai + e j aj)A k {Eiai,...,E k a k ){l + 0{nq)~* 



l<i<j<k 



For the purpose of considering averages of even forms or odd forms separately, we note 



that 
(4.5.9) 



E £ / L /(§ + a i) ■ ■ ■ L A\ + a k) 

f£H n (q) 

= X(\ + a 1 ) E h L f (\- ai )L f (\ + a 2 )...L f (\ + a k ) 

f&H 2 (q) 

= X{\ + ai)M{\- -a ls a 2 , • • • , « fc )(l + 0(nq)-^ +£ ). 
By looking at the combinations 

(4.5.10) E" L f {\ + a x )...L f {\ + a k )± E* e f L f {\ + a x ) . . . L f {\ + a k ) 

feH n (q) feH n (q) 

this leads to 

Conjecture 4.5.2. With A k (ai, . . . ,a k ) as in Coniecture \1.5.5l we have 



L f (± + a 1 )...L f (± + a k ) 



f<£H n (q) 
f even 



k k 



(4.5.11) 
and 



e j= ±lj = l 

x I I + ti(Xi + ^j^j)A k {tiai, . . . ,e k a k )(l + 0(nq)~* +s ), 

l<i<j<k 



E L f {\ + ai )...L f {\ + a k ) 



f&H n (q) 
f odd 



(4.5.12) 



ej =±ij=i 

x ] [ C(l + + e i a i )A jfc (eia;i, . . . ,e k a k )(l + 0{nq)~^ +£ ). 



3=1 



l<i<j<k 
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The above formulae can be written as contour integrals using Lemma l2.5.2[ giving ex- 
pressions analogous to those in Conjecture 11.5.51 In particular, expressing (|4.5.11j) as a 
contour integral and then letting aj — * gives Conjecture 11.5.51 

5. Numerical calculations 

We compare our conjectures with some numerical calculations. The agreement is very 
good. These calculations involve numerically approximating the coefficients in the conjec- 
tured formulae, and numerically evaluating the mean value. Both of those calculations are 
interesting and we will give more details and examples in a subsequent paper. 

5.1. Unitary: Riemann zeta-function. The coefficients of P2(x) in Conjecture 11.5. II can 
be written explicitly in terms of known constants. When k = 2 the function G{a\ 1 «2, «3, (Xj) 
that appears in Conjecture 11.5. II equals 

2 

(5.1.1) C( 2 + oli + «2 - « 3 - uaT 1 \\ C(l + Qfj - atk+j), 

which is given in Section l2~2l 
But 

(5.1.2) C(l + s) = s- 1 + 7 - 7i* + - |p + • • • 
and 



(5.1.3) C (2 + .)-' = 4 - ^ + t + 

7T^ 7T 4 7T b 



In the latter, the terms up to s 4 were evaluated using Maple. Computing the residue in 
Conjecture 11.5.11 gives 



P*(x) = ^ + (V - 3C'(2)) x 3 

+ JL (-487C'(2)7r 2 - 12C"(2)7r 2 + 7 7 2 7r 4 + 144C'(2) 2 - 2 7l 7r 4 ) x 2 

7T b 

+ ^ f 6 7 3 7r 6 - 84 7 2 C / (2)tt 4 + 24 7i C / (2)tt 4 - 1728C'(2) 3 + 5767C'(2)V 



7T° 



+ 288C'(2)C // (2)vr 2 - 8C' // (2)vr 4 - lO^vr 6 - 72 7r 6 - 48 7 C"(2)7r 4 )x 



+ ^-12C /,,/ (2)tt 6 + 3672C'(2)vr 6 + 9 7 4 vr 8 + 21 7l 2 7r 8 + 432C"(2) 2 7r 4 

+ 3456 7 C , (2)C // (2)7r 4 + 30247 2 C'(2) 2 7r 4 - 367 2 7l 7r 8 - 252 7 2 C"(2)7r 6 

+ 37727T 8 + 727iC"(2)7r 6 + 3607i7C'(2)7r 6 - 2167 3 C'(2)tt 6 

- 864 7i C / (2) 2 tt 4 + 5 73 7r 8 + 576C , (2)C / "(2)vr 4 - 20736 7 C'(2)V 

(5.1.4) - 15552C ,/ (2)C'(2)V - 967C' // (2)tt 6 + 62208C'(2) 4 
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in agreement with a result implied in the work of Heath-Brown [H-BQ (see [d where, using 
[H-BQ, the same polynomial is worked out, although there are some slight errors). Numeri- 
cally, 

p 2 (x) = 0.0506605918211688857219397316048638 x 4 

+ 0.69886988487897996984709628427658502 x 3 
+ 2.425962198846682004756575310160663 a; 2 
+ 3.227907964901254764380689851274668 x 

(5.1.5) + 1.312424385961669226168440066229978 

There are several ways that one can numerically compute the coefficients of Ps(x), and 
these will be described in a future paper. We found 

P 3 (x) = 0.000005708527034652788398376841445252313 x 9 
+ 0.00040502133088411440331215332025984 a; 8 
+ 0.011072455215246998350410400826667a; 7 
+ 0.14840073080150272680851401518774 x 6 
+ 1.0459251779054883439385323798059 a; 5 
+ 3.984385094823534724747964073429 x 4 

+ 8.60731914578120675614834763629 a; 3 + 10.274330830703446134183009522 a; 2 

(5.1.6) + 6.59391302064975810465713392 x + 0.9165155076378930590178543. 

One notices that the leading coefficient is much smaller than the lower order coefficients, 
which means that, in numerical calculations, the lower order terms will contribute signifi- 
cantly. One might suppose that the coefficients of Pk{x) are always positive. Unfortunately, 
while this is true for P\, . . . , P±, by k = 5, negative coefficients begin to appear (see Ta- 
ble ettt]i . 



Table 15.1.21 depicts 



(5.1.7) / \C(l + zt)\ 6 dt 

Jc 

as compared to 

(5.1.8) / F 3 (]og(t/27r))dfc, 

Jc 

along with their ratio, for various blocks [C, D] of length 50000, as well as a larger block of 
length 2,350,000. The data agree with our conjecture and are consistent with a remainder 
of size \D-C\^D e . 

One can also look at smoothed moments, for example, 

/>oo 

(5.1.9) / \((l + it)\ 2k exp(-t/T)dt 

Jo 

as compared to 

/■oo 

(5.1.10) / P fe (log(t/27r))exp(-t/T)di. 

Jo 
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Table EH31 compares these with T = 10000, for k = 4,3, 2, 1. 

For k = 3, 4 the data agrees to roughly half the decimal places. This supports our 
conjecture that the error term in the conjectured mean values is 0(T5 +e ). For k = 1 the 
numerics suggest Corollary 11.6.31 

5.2. Symplectic: quadratic Dirichlet L-functions. We have computed the polynomials 
Qk of Conjecture 11.5.31 for k = 1, 2, . . . , 8, separately for d < and d > 0. Table 15.2.11 lists 
these polynomials for d < 0, while in Table lB~2.2l we consider d > 0. Again notice the small 
size of the leading coefficients. 

Table IB~2.3I compares, for d < 0, conjectured moments for k = 1, . . . , 8 against numeri- 
cally computed moments, 

(5-2.1) L (hxd) k g(\d\) 

d<0 

versus 

(5.2.2) £f Q fc (log|d|)0(|d|) 

d<0 

where g is the smooth test function 

1, if < t < 850000; 



(5.2.3) g(t) 



exp I 1 - _ (t _ 8 1 50000)2 J , if 850000 < t < 1000000; 

V 1 (150000) 2 / 

0, if 1000000 < t. 



Table I5~2~4l compares the same quantities, but for d > 0. 

Figure □ depicts, for k = 1, . . . , 8 and X = 10000, 20000, . . . , 10 7 , 

(5-2.4) 

-X<d<0 

divided by 

(5.2.5) Q^°M)- 

-X<d<0 

One sees the graphs fluctuating above and below one. Interestingly, the graphs have a similar 
shape as k varies. This is explained by the fact that large values of L(l/2, \d) tend to skew 
the moments, and this gets amplified as k increases. 

Figure 121 depicts the same but for < d < X. 

The values L(|, Xd) were computed using a smoothed form of the approximate functional 
equation which expresses the L-function in terms of the incomplete Gamma-function (see 
for example |L] ). 

5.3. Orthogonal: twists of a GL(2) L-function. Let 

oo 

(5.3.i) M*) = £-fe 



n=l 
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be the L-function of conductor 11 of the elliptic curve 

(5.3.2) y 2 + y = x 3 - x 2 . 

The coefficients a n are obtained from the cusp form of weight two and level 11 given by 

oo oo 

(5.3.3) a ^ = ? IR 1 - ^ - <7 Un ) 2 - 

n=l n=l 

Expanding the right side using Euler's pentagonal theorem provides an efficient means to 
compute the a n 's. 

L n (s) satisfies an even functional equation (i.e. e = +1): 

(5.3.4) (^j T(s + l)L u (s) = (^j r(§ - s)L u (l - s), 
and may be written product over primes 

(5.3.5) Ln(s) = r TT -. 

1 - ll"*"* p % 1 - a pP - s —2 + 

Consider now quadratic twists of L u (s), 

oo 

(5.3.6) L u (s, Xd) = ^ ~rr;Xd(n). 



n=l 712 



with (d, 11) = 1. Lu(s,Xd) satisfies the functional equation 

(5.3.7) L n (s, Xd ) = X,(-H)^4 (^X ' \d\ 2 ^~ s) L U {1 - s, Xd ). 

1 {s+ 2) V 11 2 / 

We wish to look at moments of L\\{\,Xd) but only for those L(s,Xd) that have an even 
functional equation, i.e. Xd{— H) = 1- We further only look at d < since in that case a 
theorem of Kohnen and Zagier |KZj enables us to easily gather numerical data for Ln(|, Xd) 
with which to check our conjecture. 

When d < 0, Xd{— 1) = — 1, hence, in order to have even functional equation, we require 
Xd(H) = — 1) i- e - d — 2, 6, 7, 8, 10 mod 11. Conjectured formula (|4.4.22j) combined with 
Lemma f2.5.2l gives an estimate for the sum over fundamental discriminants 

(5.3.8) L n (^Xd) k = ^ T fe (log |d|) + 0(Z^ +£ ) 

-D<d<0 -D<d<0 
d=2,6,7,8,10 mod 11 d=2,6,7,8,10 mod 11 

where, as in Section [Ol is the polynomial of degree k(k — l)/2 given by the /c-fold residue 

(5 'r t) = '-"f^r^ 1-4 ^1->*D\*T.U , ^ . . . *>, 



fc! (2vrz) 
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where 

k 



(5.3.10) R xl {z u •••,*») = M*i, ••-,**) f[ ( If. * Z jl (-rsT) 2 II C(l + Zi + Zj ), 



j=i v v j/ \ / / i<i<j<k 
and A fc is the Euler product which is absolutely convergent for ^ \zj\ < \ 



2' 



(5.3.11) A h (z 1 ,...,z k )=Y[R UtP (z 1 ,...,z k ) Yl f 1 - i+! 1+ ,J 



l<i<j<fc 



with, for p 7^ 11, 
(5.3.12) 



ftn,p — I I - 
and 



1 \ / 1 1 TT 1 T-r 1 



P J y p ^ 2 y ^ 1 _ ap jr> 1 Zj • _)_ p 1 22j TT I _|_ 1 Zj _|_ p 



(5.3.13) i2n,n=n 



1 + 

i=i 

Numerically, it is more challenging to compute the polynomials TV First, using 

(5.3.14) J] C0- + Zi + Zj) 

l<i<j<k 

to estimate the sum over primes of ()4.4.17|) makes a poor approximation and one would do 
better to use the Rankin-Selberg convolution L-function of Ln(s) with itself. However, it is 
simpler to work with £, and we thus computed the first 4 moment polynomials of £ii(§, Xd) 
but to low accuracy. The coefficients of these polynomials are given to 2-5 decimal place 
accuracy in Table lo~3~Tl 

In Table l5~3.2l we compare moments computed numerically with moments estimated by 
our conjecture. The two agree to within the accuracy we have for the moment polynomial 
coefficients. We believe that if one were to compute the coefficients to higher accuracy, one 
would see an even better agreement with the data. 

While one can compute i<Li(§, Xd) using standard techniques (see jCohj ). in our case we 
can exploit a theorem of Kohnen and Zagier |KZj which relates £n(|, Xd), for fundamental 
discriminants d < 0, d = 2, 6, 7, 8, 10 mod 11, to the coefficients cn(|d|) of a weight 3/2 
modular form 

(5.3.15) L u (lxd) = K llCll (\d\) 2 /Vd 

where Ku is a constant. The weight 3/2 form in question was determined by Rodriguez- 
Villegas (private communication) 

oo 

J2 c n(n)q n = (G 1 (q) - 2 (q))/2 

n=l 

(5.3.16) = -q 3 + q i + q 11 + q 12 -q 15 -2q le -q 20 ... 
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where 

(5.3.17) 6 1 (q) = q x2+lly2+llz2 = 1 + 2g 4 + 2q u + 4g 12 + 4g 15 + 2g 16 + 4g 20 . . . 



x = y mod 2 



and 

(5.3.18) 6 2 (q) = q {x2+Uy2+33z2)/3 = 1 + 2q 3 + 2q 12 + 6g 15 + 6g 16 + 6g 20 .... 

{x ,y ,z) £.7? 
x = y mod 3 
y = z mod 2 

This was used to compute the cn(|d|)'s for d < 85, 000, 000. 

Evaluating the left side of (|5.3.8|) in a more traditional manner for d = —3, and comparing 
with the right side, we determined 

(5.3.19) K n = 2.917633233876991. 
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Table 5.1.1. Coefficients of P k {x) = c (k)x k2 + c 1 (k)x k2 1 + • • • + c k 2(k), 
for k = 4,5,6,7. Notice the relatively small size of cq(/c). We believe the 



coefficients to be correct to the number of places listed, except in the cases 
indicated by question marks, where the numerics have not quite stabilized. 
Two different methods were used to compute the coefficients. The former, for 
< r < 7, gave us higher precision but was less efficient, while the latter for 
r < 49, was more efficient but required using less precision. 
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[C, D] 


conjecture (|5. 1.8)1 


realitv (15.1.71) 


ratio 


[0,50000] 


7236872972.7 


7231005642.3 


.999189 


[50000,100000] 


15696470555.3 


15723919113.6 


1.001749 


[100000,150000] 


21568672884.1 


21536840937.9 


.998524 
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26381397608.2 


26246250354.1 


.994877 


[200000,250000] 


30556177136.5 


30692229217.8 


1.004453 


[250000,300000] 


34290291841.0 


34414329738.9 
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[300000,350000] 
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37683495193.0 
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40843941365.7 


40566252008.5 
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62765220708.6 


1.001663 
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64636649728.0 


64227164326.1 
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[900000,950000] 


66552376294.2 


65994874052.2 


.991623 


[950000,1000000] 


68412937271.4 


68961125079.8 


1.008013 


[1000000,1050000] 


70222493232.7 


70233393177.0 


1.000155 


[1050000,1100000] 


71984709805.4 


72919426905.7 


1.012985 


[1100000,1150000] 


73702836332.4 


72567024812.4 


.984589 


[1150000,1200000] 


75379769148.4 


76267763314.7 


1.011780 


[1200000,1250000] 


77018102997.5 


76750297112.6 


.996523 


[1250000,1300000] 


78620173202.6 


78315210623.9 


.996121 


[1300000,1350000] 


80188090542.5 


80320710380.9 


1.001654 


[1350000,1400000] 


81723770322.2 


80767881132.6 


.988303 


[1400000,1450000] 


83228956776.3 


83782957374.3 


1.006656 


[0,2350000] 


3317437762612.4 


3317496016044.9 


1.000017 



col- 



umn, i.e. integrals involving (, were computed using Mathematica. 



k 


(|5.1.9|) 




difference 


relative 
difference 


1 


79499.9312635 


79496.7897047 


3.14156 


3.952 x 10" 5 


2 


5088332.55512 


5088336.43654 


-3.8814 


-7.628 x 10- 7 


3 


708967359.4 


708965694.5 


1664.9 


2.348 x 10" 6 


4 


143638308513.0 


143628911646.6 


9396866.4 


6.542 x 10" 5 



C versus Conjecture 11.5.11 " 



Table 5.1.3. Smoothed moment o 
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r 


dr(l) 


d r (2) 


d r (3) 


d r (A) 





.3522211004995828 


.1238375103096e-l 


.1528376099282e-4 


.31582683324433e-9 


1 


.61755003361406 


.18074683511868 


.89682763979959e-3 


.50622013406082e-7 


2 




.3658991414081 


.17014201759477e-l 


.32520704779144e-5 


3 




-.13989539029 


.10932818306819 


.10650782552992e-3 


4 






.13585569409025 


.18657913487212e-2 


5 






-.23295091113684 


.16586741288851e-l 


6 






.47353038377966 


.59859999105052e-l 


7 

8 








.52311798496e-2 
-.1097356195 


9 








.55812532 


10 








.19185945 


r 


d r (5) 


dr(6) 


dr(7) 


d r (8) 





.671251761107e-16 


.1036004645427e-24 


.886492719e-36 


.337201e-49 


1 


.23412332535824e-13 


.67968140667178e-22 


.98944375081241e-33 


.59511917e-46 


2 


.35711692341033e-ll 


.20378083365099e-19 


.51762930260135e-30 


.500204322e-43 


3 


.31271184907852e-9 


.36980514080794e-17 


.16867245856115e-27 


.2664702284e-40 


4 


.17346173129392e-7 


.45348387982697e-15 


.38372675160809e-25 


.1010164552e-37 


5 


.63429411057027e-6 


.39728668850800e-13 


.64746354773372e-23 


.29004988867e-35 


6 


.15410644373832e-4 


.2563279107877e-ll 


.84021141030379e-21 


.65555882460e-33 


7 


.2441498848698e-3 


.12372292296e-9 


.85817644593981e-19 


.11966099802e-30 


8 


.2390928284571e-2 


.44915158297e-8 


.70024645896e-17 


.17958286298e-28 


9 


.127561073626e-l 


.1222154548e-6 


.4607034349989e-15 


.22443685425e-26 


10 


.24303820161e-l 


.2461203700e-5 


.2455973970377e-13 


.2357312577e-24 


11 


-.333141763e-l 


.3579140509e-4 


.106223013225e-ll 


.20942850060e-22 


12 


.25775611e-l 


.3597968761e-3 


.3719625461492e-10 


.15805997923e-20 


13 


.531596583 


.230207769e-2 


.1048661496741e-8 


.10159435845e-18 


14 


-.325832 


.7699469185e-2 


.2357398870407e-7 


.55665248752e-17 


15 


-1.34187 


.4281359929e-2 


.416315210727e-6 


.25985097519e-15 


16 




-.2312387714e-l 


.564739434674e-5 


.103134457e-13 


17 




.109503 


.56831273239e-4 


.346778002e-12 


18 




.2900464 


.40016131254e-3 


.982481680e-ll 


19 




-.9016 


.1755324808e-2 


.232784142e-9 


20 




-.89361 


.340409901e-2 


.456549799e-8 


21 




-.181 


-.2741804e-2 


.7309216472e-7 


22 






.353555e-3 


.9368893764e-6 


23 






.117734 


.9348804928e-5 


24 






.20714e-l 


.69517414e-4 


25 






-.9671 


.356576507e-3 


26 






-.284 


.1059852e-2 


27 






1.3 


.8242527e-3 


28 






-1. 


-.206921e-2 


29 








.181031e-l 


30 








.862815e-l 


31 








-.14025 


32 








-.91619 


33 








-.942 


34 








-.153e-l 


35 








-.3 ? 


36 








? 



Table 5.2.1. Coefficients of Q k (x) = d (k)x k ^ k+1 ^ 2 + d^x^ 1 ^ 2 + • • • , 



for k — 1 8. odd twists, d, < 0. 
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r 


e r (l) 


e r (2) 


e r (3) 


e r (4) 





.3522211004995828 


.1238375103096e-l 


.1528376099282e-4 


.31582683324433e-9 


1 


-.4889851881547 


.6403273133043e-l 


.60873553227400e-3 


.40700020814812e-7 


2 




-.403098546303 


.51895362572218e-2 


.19610356347280e-5 


3 




.878472325297 


-.20704166961612e-l 


.4187933734219e-4 


4 






-.4836560144296e-l 


.32338329823195e-3 


5 






.6305676273171 


-.7264209058150e-3 


6 






-1.23114954368 


-.97413031149e-2 


7 








.6254058547e-l 


8 








.533803934e-l 


9 








-1.125788 


10 








2.125417 


r 


e r (5) 


e r (6) 


e r (7) 


er(8) 





.671251761107e-16 


.1036004645427e-24 


.886492719e-36 


.337201e-49 


1 


.2024913313373e-13 


.6113326104277e-22 


.91146378e-33 


.556982629e-46 


2 


.261100345555e-ll 


.16322243213252e-19 


.437008961e-30 


.43686422e-43 


3 


.187088892376e-9 


.2605311255687e-17 


.1297363095e-27 


.216465856e-40 


4 


.8086250862418e-8 


.2766415183453e-15 


.2670392090e-25 


.7604817313e-38 


5 


.2126496335545e-6 


.2056437432502e-13 


.404346681e-23 


.201532781e-35 


6 


.319415704903e-5 


.10957094998959e-ll 


.46631481394e-21 


.418459324e-33 


7 


.21201987479e-4 


.42061728711797e-10 


.41831543311e-19 


.698046515e-31 


8 


-.33900555230e-4 


.11491097182922e-8 


.29548572643e-17 


.951665168e-29 


9 


-.775061385e-3 


.21545094604323e-7 


.1652770327e-15 


.1073015400e-26 


10 


.333997849e-2 


.25433712247032e-6 


.73192383650e-14 


.1008662234e-24 


11 


.22204682e-l 


.1448397731463e-5 


.25506469557e-12 


.7945270901e-23 


12 


-.1538433 


-.2179868777201e-5 


.6901276286e-ll 


.5257922143e-21 


13 


-.19794e-l 


-.54298634893e-4 


.141485467e-9 


.2924082555e-19 


14 


2.01541 


.1698771341e-3 


.210241720e-8 


.1363867915e-17 


15 


-4.451 


.22887524e-2 


.20651382e-7 


.5311448709e-16 


16 




-.1042e-l 


.101650951e-6 


.1714154659e-14 


17 




-.4339429e-l 


-.16979129e-6 


.453180963e-13 


18 




.343054 


-.37367e-5 


.9644403068e-12 


19 




-.1947171 


.97069e-5 


.160742335e-10 


20 




-3.16910 


.18351e-3 


.200188929e-9 


21 




7.31266 


-.54878e-3 


.16931900e-8 


22 






-.5621e-2 


.7257434e-8 


23 






.284e-l 


-.14329111e-7 


24 






.639e-l 


-.25913136e-6 


25 






-.7 


.6473933e-6 


26 






.86 


.138673e-4 


27 






5. 


-.2339e-4 


28 






-.Ie2 


-.48124e-3 


29 








.162e-2 


30 








.976e-2 


31 








-.83e-l 


32 








-.62e-l 


33 








2. 


34 








-2. 


35 








-9. 


36 








30.? 



Table 5.2.2. Coefficients of Q k {x) = e (k)x k( - k+i y 2 + ei {k)x k ( k+1 V 2 + • • • , 



for k = 1 8. even twists, d > D. 



INTEGRAL MOMENTS OF L-FUNCTIONS 



(>!) 



k 


reality (15.2.111 


conjecture (|5.2.2|) 


ratio 


1 


1460861.8 


1460891. 


0.99998 


2 


17225813.8 


17226897.5 


0.999937 


3 


316065502.1 


316107868.6 


0.999866 


4 


7378585496. 


7380357447.1 


0.99976 


5 


198754711593.6 


198809762196.4 


0.999723 


6 


5876732216291.7 


5877354317291.3 


0.999894 


7 


185524225881950. 


185451557119001. 


1.000392 


8 


6149876164696600 


6141908614344770 


1.0013 



Table 5.2.3. Smoothed moment of L(|,Xcz) versus Conjecture 11.5.31 for 
fundamental discriminants —1000000 < d < 0, and k = 1, . . . ,8. 



k 


realitv (|5T27T]> 


conjecture (|5.2.2|) 


ratio 


1 


1144563.5 


1144535.5 


1.000024 


2 


9252479.6 


9252229.9 


1.000027 


3 


109917867.0 


109917367.9 


1.0000045 


4 


1622521963.4 


1622508843.4 


1.0000081 


5 


27321430060. 


27320230686. 


1.000043 


6 


501621762060.6 


501542204848.7 


1.000159 


7 


9787833470714.1 


9783848274459.6 


1.000407 


8 


199831160877919 


199664775232854 


1.000833 



Table 5.2.4. Smoothed moment of versus Conjecture 11.5.31 for 

fundamental discriminants < d < 1000000, and k — 1, ... ,8. 



r 


/r(l) 


/r(2) 


/r(3) 


/r(4) 





1.2353 


.3834 


.00804 


.0000058 


1 




1.850 


.209 


.000444 


2 






1.57 


.0132 


3 






2.85 


.1919 


4 








1.381 


5 








4.41 


6 








4.3 



Table 5.3.1. Coefficients of T k (x) = / (£;)x fe(fe - 1)/2 + A(A0z fc(fe ~ 1)/2-1 + • • ., 
for k=l,2,3,4. 



k 


left side (15.3.811 


right side (I5.3.8|l 


ratio 


1 


14628043.5 


14628305. 


0.99998 


2 


100242348.8 


100263216. 


0.9998 


3 


1584067116.8 


1587623419. 


0.998 


4 


41674900434.9 


41989559937. 


0.993 



Table 5.3.2. Moments of Xd) versus their conjectured values, for fun- 

damental discriminants —85,000,000 < d < 0, d — 2, 6, 7, 8, 10 mod 11, and 
k = 1, ... ,4. The data agree with our conjectures to the accuracy to which 
we have computed the moment polynomials T^. 
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0.9998 -4«% + +^ 
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Figure 1. Horizontal axis in each graph is X. These graphs depict the first 
eight moments, sharp cutoff, of L(l/2,Xd), ~X < d < divided by the 
conjectured value, sampled at X = 10000, 20000, . . . , 10 7 . One sees the graphs 
fluctuating above and below one. Notice that the vertical scale varies from 



INTEGRAL MOMENTS OF L-FUNCTIONS 
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Figure 2. Same as the previous fingure, but for < d < X. 



